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Report  on  the  project  JV®  2088p 


“THEORETICAL  AND  EXPERIMENTAL  RESEARCH  OF  CAPABILITIES  OF  MHD 
TECHNOLOGY  TO  CONTROL  GAS  FLOW  WITH  NON-EQUILIBRIUM  IONIZATION” 

1.  INTRODUCTION 

The  purpose  of  the  project  is  theoretical  and  experimental  research  of  capabilities  of 
using  of  MHD  technology  to  control  gas  flow  with  non-equilibrium  ionization.  This  report  is 
devoted  to  theoretical  investigation  of  the  problem. 

“Cold”  gas  flows  are  considered,  such  that  thermal  ionization  of  gas  flow  is 
negligible.  The  analysis  of  various  methods  of  ionization  is  carried  out,  and  mathematical 
models  of  an  ionizer  and  MHD  generator  are  developed.  Requirements  to  ionizer,  MHD 
generator  and  flow  parameters  at  which  self-sustained  operational  mode  of  ionizer  and  MHD 
generator  is  realized  are  formulated.  Possibilities  of  using  of  MHD  control  in  gas-dynamical 
systems  are  considered. 

Traditional  using  of  MHD  transformation  of  energy  to  produce  electric  power  is  well 
understood  both  theoretically  and  experimentally  nowadays.  Nontraditional  using  of  MHD  to 
control  flow  parameters  and  hence  performance  of  gas-dynamical  systems,  which  include 
MHD  system  as  a  subsystem,  is  not  well  understood.  MHD  interaction  acts  as  bulk  force  and 
power  upon  flow  and  allows  varying  flow  parameters.  In  this  case  no  mechanical  parts  are 
used  to  control  the  flow  parameters.  So  this  method  of  control  can  be  considered  as  a 
perspective  one  for  gas-dynamical  system  with  a  fixed  geometry  especially  for  high  enthalpy 
flows.  The  main  requirements  to  realize  MHD  influence  on  flow  are  the  flow  conductivity 
ensuring.  In  the  most  of  gas-dynamical  systems  thermal  ionization  of  flow  is  inappreciable 
and  therefore  it  is  necessary  to  use  methods  of  non-equilibrium  plasma  creation.  To  ensure 
autonomous  regime  of  system  of  MHD  control  self- sustained  operational  mode  of  ionizer  and 
MHD  generator  is  considered.  For  the  operational  mode  power  spent  on  flow  ionization  is 
less  than  power  produced  by  MHD  generator. 

As  a  possible  range  of  application  for  MHD  control,  two  types  of  gas-dynamical 
system  can  be  considered:  a)  gas-dynamical  systems  which  must  ensure  variable  parameters 
of  flow  at  outlet  of  system  with  fixed  parameters  at  inlet  of  one  and  b)  gas-dynamical  systems 
which  must  ensure  limited  range  of  variation  for  some  flow  parameters  at  outlet  of  system 
with  variable  parameters  at  inlet. 

In  the  project  we  consider  application  of  MHD  systems  for  scramjet  control. 

2.  MODEL  OF  MHD  GENERATOR  WITH  NONEQUILIBRIUM  CONDUCTIVITY. 

One  of  the  potential  applications  of  MHD  generator  with  non-equilibrium 
conductivity  is  flow  control  in  scramjet,  which  allows  us  to  improve  scramjet  performance. 
Scramjet  with  MHD  control  under  “AJAX”  concept  use  MHD  generators  that  are  located 
upstream  the  combustion  chamber.  Thus  static  temperature  of  gas  in  the  MHD  generators 
must  be  not  greater  than  two  thousand  degrees  of  Kelvin.  Equilibrium  ionization  of  gas  at 
these  conditions  is  negligible,  so  to  provide  noticeable  MHD  effect  we  must  use  methods  of 
nonequilibrium  plasma  creation.  For  this  purpose  it  is  possible  to  use  different  types  of  gas 
discharge,  microwave  radiation  or  beams  of  fast  charged  particles.  In  all  the  methods  we  need 
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In  papers  [8,  9]  analytical  solution  for  set  of  equations  (1)  for  Faraday  MHD  generator 
was  obtained  in  the  form  convenient  for  analysis  of  complex  gas-dynamic  systems  including 
MHD  systems.  In  [10]  results  of  [8,  9]  were  extended  on  the  case  of  Hall  MHD  generator.  To 
develop  mathematical  model  of  MHD  generator  with  nonequilibrium  conductivity  we  use 
results  from  papers  [8-10]  taking  into  account  model  of  nonequilibrium  plasma  developed 
above. 


According  to  [10]  relations  of  flow  parameters  at  MHD  generator  exit  (subindex  2)  to 
the  flow  parameters  at  MHD  generator  entrance  (subindex  1)  can  be  written  in  the  form: 
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where  £  is  factor  which  is  determined  by  the  flow  regime,  7  is  the  enthalpy  extraction  ratio, 
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,  Mi  is  the  Mach  number  at  the  MHD  generator  entrance,  y  is  the 
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Table  3 

Values  of  £  parameter  for  particular  regimes  of  MHD  flow 


Flow 

regime 

Type  of  MHD  generator 

Faraday 

Hall 

/>const 

y-1 

y-1 

p- const 

0 

0 

7-const 

-1 

-1 

M-const 

1  k 
(l-k)co 

!  Pk(\-k) 

(l  +  p2k2)(D 

v-const 

1 

1  -k 

1  +  P2k 

1  +  P2k2 

Values  of  parameter  £  for  some  particular  flow  regimes  are  presented  in  the  Table  3. 
Value  <;=0  corresponds  to  flow  with  constant  static  pressure  along  the  channel  length.  When 
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p> 0 ,  static  pressure  increases  along  the  channel  length.  Cross-sectional  area  of  MHD  channel 
A  as  a  function  of  coordinate  x  can  be  determined  from  condition  /?vd=const.  Function  of  A 
from  the  enthalpy  extraction  ratio  at  given  value  of  factor  £  can  be  easily  obtained  by  using 


formulae  (17): 


A  =  AJ  G(ip,  4,  Mx ,  ? ])DJZ(y/,  J,  Mx ,  77) 


.  According  to  [10],  the  case  of 


constant  cross-sectional  area  of  the  MHD  generator  d=const  is  described  well  when 


£  =  (ys) .  The  (if/')  function  is  determined  by  the  expression:  ^{y)  =  y- 1  + 


r+w 

M?  -1 


.  This 


ratio  is  obtained  from  condition  dAI  drj  =  0 . 

To  obtain  dependence  of  the  enthalpy  extraction  ratio  upon  the  length  of  MHD 
generator  for  a  given  value  of  d,  parameter  it  is  necessary  to  use  the  next  relation  from  paper 
[8]: 

^  =  ± (19) 
dT  dT  %(\-kfcrB2v 

Nonequilibrium  conductivity,  according  to  (7-10),  depends  upon  flow  parameters  and 
power  density  spent  on  flow  ionization  in  the  form: 
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By  substituting  (17,20)  in  (19)  we  obtain  the  next  expression  for  MHD  interaction  parameter 
Sv  as  a  function  of  the  enthalpy  extraction  ratio 
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X  = 


1  for  Faraday  MHD  generator 

fi2 1 (l  +  /?2  j  for  Hall  MHD  generator 


where  Sv  =  a^B2  Lj p^\\  .  It  is  supposed  that  power  density  spent  on  flow  ionization  is  constant 
in  channel  of  MHD  generator  (qi=q  1),  cr,  =  a(qt,  .  Another  physically  reasonable  situation 


can  be  determined  by  a  constant  value  of  e-beam  current  density  along  the  channel  length.  In 
this  case  the  power  density  spent  on  ionization,  according  to  (6),  is  proportional  to  gas 
density  qt~  p .  Relation  between  Sv  and  rj  in  this  approach  will  takes  up  the  form: 
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Formulae  (21)  determines  MHD  interaction  parameter  as  a  function  of  enthalpy  extraction 
ratio  Sv(r/)  in  approaches  with  a  constant  power  density  spent  on  ionization  (21a)  and  with  a 
constant  e-beam  current  density  (21b).  By  inverting  the  relations  we  obtain  dependence  of 

V(Sv). 

So,  set  of  formulae  (17,18,21)  allows  us  to  determine  flow  parameters  along  the 
length  of  MHD  generator  with  non-equilibrium  conductivity.  In  common  case,  according  the 
formulae  (7-1 1),  the  flow  parameters  change  causes  the  critical  parameters  to  be  modified.  In 
order  to  take  into  account  these  non-local  effects  on  the  limits  of  the  self-sustained 
operational  mode  for  MHD  generator  with  nonequilibrium  conductivity  we  introduce  the  next 
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Fig.3.  Electron  concentration  in  air  as  a  function  of  power  loss  by  e-beam  in  unit  of  volume. 


Fig.4.  Ionization  fraction  in  air  as  a  function  of  power  loss  by  e-beam  in  unit  of  volume. 
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Fig.5  Ionization  fraction  in  air  plasma  as  a  function  of  power  loss  by  e-beam  in  unit  of  volume. 


Fig.6.  Power  density  produced  by  MHD  generator  except  for  power  density  input  to  ionization 
(qg-qt)  as  a  function  of  qt  and  B.  MHD  generator  is  located  in  cross-section  2  (Fig.l),  Moo=6. 
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Fig.9  Relative  value  of  power  put  to  flow  ionization  in  channel  of  MHD  generator  as  a 
function  of  the  enthalpy  extraction  ratio  for  various  values  of  Mach  number  and  load  factor. 
(qi  -  constant) 

a)  Faraday  MHD  generator;  b)  Hall  MHD  generator  with  J3=  2 
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Fig.  10  Relative  power  characteristics  of  MHD  generator  with  nonequilibrium  conductivity  as  a 
function  of  MHD  interaction  parameter.  Solid  curves  for  Faraday  MHD  generator  and  dashed 
curves  for  Hall  MHD  generator  with  p=2.  Mf=  3.3,  {ql  -  constant)  q/qcr=0.05,  2=0.5,  /c=0.5. 
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Fig.  11  Relative  power  excess  produced  by  Faraday  MHD  generator  with  nonequilibrium 
conductivity:  k= 0.5,  Mj=  4.2  for  solid  lines  Mj= 3.3  for  dashed  lines;  black  lines  correspond  to 
constant  power  density  spent  on  ionization,  red  lines  correspond  to  constant  e-beam  current 
density,  a)  £  =  |  (i//) ,  b)  %=() 
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Fig.  12  Relative  power  excess  produced  by  Faraday  MHD  generator  with  nonequilibrium 
conductivity:  £=0.5,  £  =  |(^);  solid  lines  M;=4.2,  dashed  lines  M;= 3.3;  black  lines 

correspond  to  constant  power  density  spent  on  ionization,  red  lines  correspond  to  constant  e- 
beam  current  density. 
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Fig.  18  Maximal  value  of  flight  Mach  number  at  which  MHD  interaction  increase  scramjet 
specific  impulse  vs.  £ parameter  for  various  values  of  y/ and  On. 


Fig.  19  Maximal  value  of  flight  Mach  number  at  which  MHD  interaction  increase  scramjet 
specific  impulse  vs.  5 parameter  for  various  y/.  For  combustion  chamber  working  in  a  mode  with 
a  constant  density  (solid  lines)  and  with  a  constant  pressure  (dashed  lines). 
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a) 


Fig.21  Relative  specific  impulse  (a)  and  relative  power  excess  (b)  vs.  MHD  interaction 
parameter  for  scramjet  with  Faraday  (solid  lines)  and  Hall  (dashed  lines)  MHD  generators. 
Red  lines  correspond  to  constant  power  density  spent  on  ionization;  black  lines  correspond  to 
constant  density  of  e-beam  current.  0n=  15°,  Mj=  10,  Fth= 0.12,  Mq= 6,  q/qCr= 0.01,  kj= 0.5 
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Fig.26  Relative  specific  impulse  for  scramjet  with  MHD  control  vs  MHD  interaction 
parameter  for  various  inlet  configurations  and  various  load  factors.  Mj=  10  F,/,=0.1 ,  Mo= 6, 
qi/qcr= 0.05,  kj=  1.2.  (a)  Faraday  MHD  generator,  (b)  Hall  MHD  generator  with  p=2 
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b) 

Fig.28  Relative  specific  impulse  (a)  and  power  excess  (b)  vs  MHD  interaction  parameter  for 
scramjet  with  Hall  MHD  generator  at  various  values  of  Hall  parameter.  Md=  10  F;/,=0.12,  Mo=6, 
q/qCr= 0.05,  ks=  1.2.  Red  curves  correspond  to  Faraday  MHD  generator. 
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Fig.30  Contours  of  relative  specific  impulse  for  scramjet  with  MHD  control.  Mo=6,  Md=  10, 
6n=\5°,  Fth=0.\2,ks=\.\,  p=2 
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a)  b) 

Fig.33  Contours  of  relative  specific  impulse  (a)  and  power  excess  (b)  for  scramjet  with  Faraday 
MHD  generator.  0N=15°,  M(i=\ 0,  F,/,=0. 1 2.  M0= 6,  L=2m,  B=5T,  ks=l.l 
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MHD  generator  with  nonequilibrium  conductivity  operates  in  self-sustained  operational 
mode  if  0<qi<qcr.  When  the  volume  V  tends  to  zero  volume,  the  value  qcr,  determined  from 
the  ratio  tends  to  value  determined  from  (12b)  or  (24)  in  local  approach.  The  calculations  of 
dependencies  of  qcr  upon  the  longitudinal  size  of  plasma  formation  were  carried  out  for  shape 
of  ionized  region  presented  in  the  Fig.37b.  We  have  considered  height  of  the  region  7?=lm. 
Position  of  the  region  is  determined  by  condition:  y<Y2  and  xi<x<X2.  Results  presented  in  the 
Fig.38  correspond  to  ionized  region  with  x/=0.  In  calculations  we  vary  length  of  the  plasma 
formation  (x^),  load  factor  k  and  angle  of  magnetic  field  orientation  a.  All  critical  powers  qcr 
are  normalized  on  the  critical  powers  obtained  for  Xi=0.5  m  and  a=90  °  for  corresponding 
values  of  the  load  factor.  Calculations  have  been  done  for  scramjet  configuration  shown  in 
the  Fig.34.  Fig.38a  demonstrates  that  critical  power  decreases  while  increasing  the  plasma 
length  practically  in  similar  way  for  various  load  factors.  Fig.38b  demonstrates  dependencies 
of  the  critical  power  upon  the  angle  a  for  identical  plasma  regions  with  X2=0.5  m.  One  can 
see  that  maximal  critical  power  is  achieved  at  a=90  °.  The  greater  is  the  load  factor  the  more 
sensitive  is  the  critical  power  value  to  the  angle  of  magnetic  field  orientation.  So  one  can  see 
that  critical  power  essentially  depends  upon  considered  configuration.  In  studying  the  MHD 
controlled  scramjet  we  verify  condition  (39)  in  all  computations.  Situations  in  which  the 
condition  is  failed  are  not  considered  in  the  report. 

The  first  step  of  our  research  is  devoted  to  investigation  of  possibilities  of  MHD 
generator  to  move  shock  wave  position.  We  consider  flow  with  initial  Mach  number  Mo=8 
and  free  stream  dynamic  pressure  40  kPa  flowing  around  the  wedge  with  the  angle  10°. 
Pressure  colour  map  without  MHD  interaction  is  shown  in  the  Fig.39a.  We  suppose  that  in 
configuring  the  magnetic  field  in  such  a  way  that  the  Lorentz  force  will  be  directed  away 
from  the  wedge  (or  will  be  positive  in  our  coordinate  system  (Fig.34))  the  angle  of  oblique 
shock  will  be  increased.  To  decrease  the  angle  of  oblique  shock  it  is  necessary  to  form  the 
Lorentz  force  directed  toward  the  wedge.  It  can  be  done  only  when  fy< 0.  We  assume  that 
By>0.  For  MHD  generator  qg< 0  and  force  fx< 0.  Thus  it  is  evidently  from  (38)  that  force  fy< 0 
only  when  Bx< 0.  So,  it  is  necessary  to  have  magnetic  field  with  a>90°  in  order  to  decrease 
the  angle  of  oblique  shock. 

Pressure  colour  map  in  the  flow  around  the  wedge  with  MHD  control  with  magnetic 
field  which  is  characterized  by  a=45°  is  shown  in  the  Fig.39b.  Range  of  MHD  interaction  is 
band  located  in  2<x<2.5  (which  is  schematically  shown  in  the  Fig  37a).  One  can  see  that 
MHD  interaction  essentially  increases  the  angle  of  the  oblique  shock.  Fig.39c  demonstrates 
pressure  colours  map  for  MHD  control  which  is  characterized  by  B=3T  and  a=90°.  Despite 
of  the  y-component  of  the  Lorentz  force  is  equal  to  zero  in  this  case,  the  oblique  shock 
deviation  is  greater  than  in  the  case  presented  in  the  Fig.39b.  In  the  Fig.  40  the  relative 
pressure  contours  (pressure  normalized  on  pressure  in  free  stream)  in  stream  flowing  around 
the  wedge  for  two  configuration  of  magnetic  field  are  compared.  One  can  see  that  in  the  case 
with  a=90°  the  MHD  interaction  provides  noticeable  increase  of  static  pressure.  In  the  Fig. 
41  the  relative  temperature  (temperature  normalized  on  temperature  in  free  stream)  contours 
in  stream  flowing  around  the  wedge  for  two  configuration  of  magnetic  field  are  compared. 
One  can  see  that  in  the  case  with  a=45°  the  MHD  interaction  practically  doesn’t  increase  the 
temperature  behind  the  shock.  So  we  can  conclude  that  increase  of  the  angle  of  oblique  shock 
in  result  of  MHD  influence  on  flow  can  be  realized  both  with  slight  modification  of  pressure 
and  temperature  behind  the  shock  (a=45°)  and  with  considerable  increase  of  these  parameters 
behind  the  shock  (a=90°).  Thus  if  we  are  interested  in  increasing  the  pressure  and 
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follows  from  the  Fig.48  that  density  distribution  tends  to  the  distribution  at  design  conditions 
while  increasing  the  power  density  spent  on  flow  ionization. 

And  now  we  will  investigate  possibilities  of  MHD  generator  to  control  flow  in 
scramjet  inlet  when  flight  Mach  number  is  less  than  design  Mach  number.  Flow  field  in 
scramjet  inlet  in  this  situation  is  shown  in  the  Fig.49.  The  serious  disadvantages  of  the  flight 
regime  are  decrease  of  flow  compression  and  air  capture  in  scramjet.  The  aim  of  our 
investigation  is  finding  the  regimes  of  MHD  influence  on  flow,  which  allow  us  to  increase 
both,  flow  compression  and  air  mass  flow  rate  in  scramjet  for  such  flight  regimes.  In 
analyzing  we  consider  two  characteristics  of  inlet:  relative  pressure  increase  (flow 
compression)  pi/po,  where  po  is  the  incoming  flow  static  pressure  and  pi  is  the  averaged  static 
pressure  at  exit  of  inlet  (in  our  computations  x=7.1 );  and  relative  air  mass  flow  (p  =  m/m0 
which  corresponds  to  ratio  of  air  mass  flow  in  inlet  m  (which  is  computed  at  exit  of  inlet 
m  =  |  pvxdy  •  Az )  to  the  free  stream  mass  flow  m0  which  is  restricted  by  the  non-disturbed 

line  of  flow  closed  in  the  cowl  lip  and  x  axis.  We  consider  only  zero  angle  of  attack,  thus 
m0  =  p0v0Yc  ■  Az  ,  here  Az  is  characteristic  width  of  inlet  in  z  direction,  Yc  is  y  coordinate  of 

the  cowl  lip  (in  our  consideration  Yc=  lm  (see  Fig.49)). 

Fig.50  demonstrates  dependencies  of  relative  mass  flow  and  flow  compression  in 
MHD  controlled  inlet  upon  power  density  spent  on  flow  ionization  for  various  configurations 
of  magnetic  fields.  Ionization  range  is  a  band  which  is  located  in  x;<x^c;+zlx.  For  Fig.50a 
x/=0  and  for  Fig.50b  x/=l ,  Ax=  lm  in  all  the  cases.  Green  lines  in  the  figure  correspond  to 
relative  mass  flow  (solid  lines)  and  relative  pressure  increase  (dashed  lines)  in  scramjet 
without  MHD  interaction.  One  can  see  that  for  x/=l  and  a=135°  there  is  a  range  of  q, 
variation  for  which  MHD  interaction  leads  to  insignificant  increase  of  air  mass  flow.  Relative 
pressure  rises  while  increasing  the  q-,  value  when  a=135°.  If  a<90°,  the  increase  of  power 
density  spent  on  flow  ionization  causes  both  the  relative  pressure  and  relative  air  mass  flow 
to  be  decreased.  This  effect  correlates  with  results  obtained  in  investigations  of  wedge  in 
magnetic  field  and  can  be  explained  by  Lorentz  force  direction. 

Fig.  51  demonstrates  dependencies  of  relative  air  mass  flow  (a)  and  relative  pressure 
increase  (b)  in  MHD  controlled  inlet  vs.  power  densities  spent  on  flow  ionization  for  various 
positions  of  ionized  region.  Orientation  of  magnetic  field  a=135°.  One  can  see  that  maximal 
air  mass  flow  can  be  achieved  when  xy=3.  In  considered  conditions  MHD  interaction 
monotonically  increase  the  pressure  in  MHD  controlled  inlet.  As  for  air  mass  flow  there  are 
some  optimal  q  -,  at  which  maximal  air  capture  in  the  inlet  can  be  realized.  Additional  increase 
of  MHD  interaction  will  provide  negative  influence  on  the  air  mass  flow. 

Fig. 5 2  demonstrates  Mach  number  contours  in  MHD  controlled  inlet  at  flight  Mach 
number  Mo= 8  which  is  less  than  design  Mach  number  M<y=10.  Flow  fields  correspond  to  three 
various  configurations  of  magnetic  field  with  magnetic  induction  5=  IT.  Level  of  MHD 
interaction  corresponding  to  these  figures  is  not  very  strong.  The  angles  of  oblique  shocks 
slightly  increase  for  a<90°.  For  a=135°  shock  wave  positions  is  practically  not  changed. 

Fig. 53  demonstrates  Mach  number  contours  in  MHD  controlled  inlet  with  magnetic 
induction  B= 3T  at  conditions  which  are  same  to  the  Fig.52.  In  this  case  the  level  of  MHD 
interaction  is  strong  enough  and  modification  of  flow  field  is  more  significant  than  in  Fig.52. 

In  order  to  view  relative  influence  of  MHD  control  on  the  air  mass  flow  and  on  the 
flow  compression  in  inlet  we  will  consider  normalized  ratios:  the  normalized  relative  air  mass 
flow  (cp/cpo),  where  cpo  is  the  relative  air  mass  flow  in  scramjet  without  MHD  control,  and 
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Fig.  3  5  Pressure  contours  in  inlet  at  design  conditions  (a)  and  fragment  of  computational  grid  (b) 
in  red  rectangular  area  marked  in  (a).  Md=  10,  Mq=  10,  6n=  15°,  F,/,=0.12 
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Fig.38  Normalized  critical  power  input  to  ionization  in  scramjet  inlet  vs  length  of  ionized 
region  (a)  and  vs  orientation  of  magnetic  field  (b).  Mo=6,  B=1T.  a=90°  in  (a)  and  x_i=0.5  m  in 
(b). 
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Fig.40  Pressure  contours  in  stream  flowing  the  wedge  with  0=10°.  Mo=8,  g;=4W/cm  ,  k= 0.5. 
a)  B=3T,  a=45°,  b)  B=3T,  a=90°.  Ionized  region  located  in  2<x<2.5. 
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Fig.54.  Normalized  relative  air  mass  flow  (cp/cpo)  and  flow  compression  (tc/tco)  in  MHD 
controlled  inlet  vs  magnetic  field  orientation  for  various  values  of  magnetic  induction. 

M0= 8,  Md=  10,  6^=15°,  Fth= 0.12,  k= 0.5,  qr\  W/cm3.  Ionized  region  is  located  in  3<x<4. 
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Fig.55.  Normalized  relative  mass-flow  rate  in  MHD  controlled  inlet  vs  load  factor  for  various 
a.  Mo= 8,  M(/=10,  <9v=15°,  F,/,=0.12,  B= 3T,  qt=  lW/cm3.  Ionized  region  is  located  in  3<x<4. 
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Fig.58.  Normalized  relative  pressure  increase  in  MHD  controlled  inlet  vs  load  factor  for 
various  a.  Mq=6,  M(/=10,  9n=  15°,  F,/,=0.12,  B= 3T,  g;=lW/cm3.  Ionized  region  is  located  in 
3<k<4. 


Fig.59  Normalized  relative  air  mass  flow  in  MHD  controlled  inlet  vs  X}  for  various  height  of 
ionized  region  Y2.  Mo= 8,  Mj=  10,  6V=15°,F(/,=0.12,  B= 3T,  a=135°,  ^,=0.1W/cm3,  k=0.5. 
Ionized  region  is  located  in  Xi<x<X/+ 1  andyXYj 
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Fig.60.  Normalized  relative  pressure  increase  in  MHD  controlled  inlet  vs  Xi  for  various 

height  of  ionized  region  Y2.  M0= 8,  Md=  10,  dN=\5°,Fth=0.\2,  B= 3T,  a=135°,  #=0.1W/ 2.0281  0  0  12.0281 
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some  energy  to  be  put  in  flow.  It  is  evidently  that  non-equilibrium  conductivity  can  be 
realized  when  energy  spent  on  flow  ionization  is  less  than  energy  produced  by  MHD 
generator.  Such  operational  mode  of  MHD  generator  and  ionizer  is  called  self-sustained 
operational  mode.  One  of  the  main  requirements  to  ionizer  is  minimal  power  spent  to 
producing  of  necessary  conductivity  of  flow.  Analysis  of  known  methods  of  ionization  shows 
that  electron  beam  is  optimal  ionizer  from  the  point  of  view  of  minimization  of  the  power 
spent  on  ionization.  According  to  [1],  electrons  with  the  energies  greater  than  1  keV  spend 
only  W,= 34  eV  to  produce  electron  and  ion  pair  in  air.  Ionization  cost  Wl  in  this  case  only  few 
times  greater  than  the  ionization  energy  of  molecules  of  air.  High-voltage  pulses  discharge, 
according  to  [2],  for  the  ratio  of  electric  field  strength  to  the  gas  number  density  E/N>5- 1 0‘15 
V-cm2,  is  characterized  by  ionization  cost  W,x66e  V.  In  principle  this  discharge  can  be 
considered  as  alternative  to  e-beam. 

Table  1 

Electron  concentration  and  conductivity  in  air  plasma  sustained  by  e-beam  with  current 
density  jb=0.1A/cm2,  (computations  of  [3]). 


e-beam  energy 
(keV) 

Pressure 

(atm) 

Temperature 

(°K) 

ne ,  10 13  cm  3 

a,  mho/m 

10 

0.013 

300 

0.68 

23.1 

600 

0.54 

33.4 

1000 

0.45 

44.4 

0.066 

300 

1.2 

10.7 

600 

1.0 

16.1 

1000 

0.91 

22.5 

0.132 

300 

1.5 

7.07 

600 

1.3 

11.6 

1000 

1.2 

16.2 

50 

0.013 

600 

0.26 

19.8 

0.066 

600 

0.50 

9.6 

0.132 

600 

0.64 

6.7 

100 

0.013 

600 

0.20 

15.9 

0.066 

600 

0.38 

7.8 

0.132 

600 

0.47 

5.3 

200 

0.013 

600 

0.16 

13.5 

0.066 

600 

0.31 

6.7 

0.132 

600 

0.38 

4.5 

To  calculate  non-equilibrium  conductivity  of  flow  in  channel  of  MHD  generator  with 
e-beam  ionizer  in  [3]  it  was  developed  model  of  ionization  of  air.  In  the  model  more  than  40 
plasma  components  and  more  than  230  reactions  of  the  plasma  components  were  taking  into 
account.  Part  of  results  computed  in  [3]  is  presented  in  the  Table  1.  MHD  generator  with  non¬ 
equilibrium  conductivity  is  a  part  of  complicated  system  -  scramjet  with  MHD  systems.  In 
order  to  analyze  the  complex  system  in  a  wide  range  of  parameters  variation,  to  formulate 
requirements  for  MHD  generator  and  ionizer,  to  determine  optimal  operational  regimes  of  the 
subsystems  we  need  to  develop  simpler  model  of  ionization  of  air.  Results  from  [3]  will  be 
used  to  obtain  simple  approximation  function  to  calculation  of  concentration  of  electrons  and 
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Concentration  of  electron  in  nonequilibrium  plasma  sustained  by  e-beam  can  be 
calculated  numerically  from  the  set  of  equations  (1).  To  estimate  energy  balance  in  MHD 
generator  with  non-equilibrium  ionization  we  consider  spatially  homogeneous  and  steady 
state  plasma.  Conductivity  of  plasma  is  determined  by  the  following  ratio 

cr  =  e2ne/(menkc),  (7) 

where  me  is  the  electron  mass,  kc  is  the  electron  scattering  rate  constant.  The  rate  constant  kc 
is  determined  by  the  ratio:  kc  =  gen  (Te)-Qen  (Te  ) ,  where  gen  is  averaged  relative  velocity  of 


electrons,  Qen  is  averaged  transport  cross-section  of  electron  collisions. 


In  the  case  when  dissociative  recombination  (2)  is  dominant  process  in  the  set  of 
equations  (1)  the  electron  concentration  ne  can  be  approximately  determined  by  the  ratio: 


n 


e 


(8a) 


Ionization  fraction  in  the  case  is  determined  by  the  ratio 

ne  „  \(khl_ 


n  \Wtp 


(8b) 


To  obtain  simple  formulas  for  calculation  of  electron  concentration  and  determine  the 
range  of  their  applicability  we  use  three  approaches: 

The  first  approach: 

Electron  concentration  is  determined  numerically  as  a  solution  of  the  set  of  equations 
(1)  with  rate  coefficients  determined  in  (2-5). 

The  second  approach: 

The  electron  concentration  will  be  determined  by  approximation  function,  which  is 
similar  to  (8a):  ne  =  at  •  ( qt  f' ,  where  coefficients  ai  and  bi  are  fitting  parameters. 


The  third  approach: 

The  ionization  fraction  will  be  determined  by  approximation  function,  which  is 

similar  to  (8b):  nel n  =  a2-(qi  / 2 ,  where  coefficients  a:  and  bj  are  fitting  parameters.  The 

fitting  parameters  a,  and  bi  for  the  second  and  the  third  approaches  are  determined  by  the 
least-squares  method  by  comparing  the  numerical  results  from  [3]  with  approximation 
functions  proposed  above. 

Electron  concentration  and  plasma  conductivity  in  MHD  generator  located  in  scramjet 
are  the  functions  of  e-beam  parameters  and  flow  parameters.  We  have  investigated  several 
configurations  for  scramjet.  Configuration  of  scramjet  for  Mach  designed  number  Mj=  10 
and  total  turning  angle  0n=  15°  is  shown  in  Fig.l.  Characteristic  cross-sections  of  scramjet  are 
denoted  in  the  figure.  Flow  parameters  are  calculated  numerically  in  2-D  Euler  approach 
(chapter  4  of  the  project).  To  provide  estimations  of  conductivity  in  channel  (cross-section  3) 
flow  parameters  are  averaged  across  the  channel.  Results  presented  in  Fig.  16  (chapter  3)  were 
used  in  calculations.  The  basic  characteristic  of  ionizer  is  power  spent  on  ionization.  So  we 
characterize  e-beam  by  integral  parameter  qt,  which  is  expressed  through  parameters  of  e- 
beam  and  flow  by  the  ratio:  qt  =  ( jje)  pY  ( Eh ) .  Fig  2  shows  dependencies  of  q,  upon  e- 


beam  current  density  in  characteristic  cross-sections  of  scramjet  for  flight  Mach  number 
Moo=  8  with  free-stream  dynamic  pressure  40kPa.  One  can  see  that  change  the  cross-section 
location  from  first  to  third  position  causes  the  ty,  value  to  be  increased. 
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Fig.3  shows  dependencies  of  electron  concentration  calculated  in  [3]  upon  g, 
magnitude.  Approximation  function  (red  line)  obtained  for  the  results  is  given  by  the  ratio: 

ne  =1.124-10127^-cm"3,  (9) 

where  qt  is  expressed  in  W/cm3.  Function  (9)  coincides  with  function  (8a)  at  corresponding 
value  of  parameter  //=  1.455 -lO^cmV1,  maximal  deviation  of  results  of  [3]  from  the 
approximation  function  is  less  than  30%.  It  follows  from  Fig.3  that  the  deviations  are 
systematic.  Thus  using  of  function  (9)  in  out  of  calculated  ranges  can  result  in  more  serious 
deviation.  Blue  lines  in  Fig.3  are  steady  state  solutions  of  (1).  In  the  calculations  we  suppose 
that  Te=Tg.  These  curves  are  in  good  agreement  with  results  of  [3],  Thus  the  second 
approach  (formula  (9))  can  be  used  to  rough  estimations  of  electron  concentration  in  the  case 
when  plasma  parameters  not  dramatically  differ  from  parameters  in  the  Fig.3.  The  first 
approach  needs  to  be  used  to  calculate  electron  concentration  more  precisely. 

Fig.4  shows  dependencies  of  ionization  fraction  of  air  plasma  upon  factor  qjn1 . 
Approximation  function  (red  line)  obtained  for  the  results  is  given  by  the  ratio: 

—  =  1.17*1 0-5  -(l034  -qi/n2)06 ,  (10) 

n  v  ' 

where  q-,  is  expressed  in  W/cm3,  n  in  cm'3.  Approximation  function  (10)  is  distinguished  from 
the  function  (8b)  by  the  index  of  power.  Function  (10)  more  precisely,  than  function  (9), 
agrees  with  numerical  results  of  [3],  Maximal  deviation  of  results  of  [3]  from  approximation 
function  (10)  is  less  than  9%.  Blue  lines  in  Fig.4  are  steady  state  solutions  of  (1).  It  follows 
from  Fig.4  that  all  dependencies  of  ionization  fraction  are  practically  coincides  in  the  range 
3-10“6  <  njn  <3-l0~5 . 

To  determine  range  of  applicability  of  formula  (10)  (the  third  approach)  to  calculation 
of  electron  concentration  in  MHD  generator  located  in  scramjet  we  have  computed  ne  in  the 
first  approach  for  all  possible  locations  of  MHD  generator  in  scramjet,  denoted  in  the  Fig.  1 , 
in  the  range  of  flight  Mach  numbers  Moo=6-hl0  at  free  stream  dynamic  pressure  40kPa.  All 
obtained  results  are  involved  in  the  area  restricted  by  the  blue  lines,  which  are  shown  in  the 
Fig.5.  Red  line  is  result  obtained  in  the  third  approach  (formula  (10)).  It  follows  from  Fig.5 
that  all  dependencies  of  ionization  fraction  are  neighbour  in  the  range  of 

0.02  <  (q.  /  n2 )  TO34  <  20  .  In  this  range  the  formula  (10)  can  be  used  with  high  reliability  to 

analyze  characteristics  of  MHD  generator  with  non-equilibrium  conductivity. 

Obtained  results  allow  us  to  estimate  power  characteristics  of  the  MHD  generator. 
Power  density  produced  by  MHD  generator  with  nonequilibrium  conductivity  qg  is 
determined  by  the  ratio: 

qg=k(l-k)cr(qi)B2v2 ,  (11) 

where  k  is  the  load  factor,  B  is  the  magnetic  induction,  v  is  the  flow  velocity.  Fig.6  shows 
dependencies  of  power  density  produced  by  MHD  generator  except  for  power  density  put 
into  flow  ionization  ( qg-qt )  as  a  function  of  qt  and  B.  One  can  see  that  there  are  limiting  values 
of  magnetic  induction  B=Bcr  for  which  qg>qj.  Increase  of  value  qt  causes  the  critical  magnetic 
induction  to  be  increased.  Increase  of  magnetic  induction  at  fixed  q,  causes  the  power 
produced  by  MHD  generator  to  be  increased.  Dependence  of  ( qg-q,)  on  the  q-,  value  is  not 
monotonic.  There  is  optimal  value  of  qh  for  which  {qg-qt)  has  maximum  at  given  value  of 
magnetic  induction. 
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To  determine  limits  at  which  self-sustained  operational  mode  of  ionizer  and  MHD 
generator  exists  we  need  to  compare  two  values:  the  power  spent  on  flow  ionisation  q„  and 
the  power  produced  by  MHD  generator  qg.  It  is  evidently  from  formulae  (7,10,1 1)  that  qg  is 
not  linear  function  of  c/(.  In  this  case  for  any  given  parameters  of  B  and  v  there  are  range  of  qt, 
which  provide  self-sustained  operational  mode  for  MHD  generator  and  ionizer.  At  fixed 
value  of  qt  we  can  determine  critical  value  (or  limiting  value)  of  magnetic  induction  Bcr  for 
which  qg=qi.  Self-sustained  operational  mode  of  MHD  generator  and  ionizer  (qg>qi)  is 
realized  when  B>Bcr.  We  have  calculated  conductivity  for  various  values  of  q,  in  approaches 
considered  above.  Function  of  critical  magnetic  induction  upon  power  spent  on  flow 
ionization  Bcr(qi)  is  determined  by  solving  the  evident  equation  qt  =  k(\-k)<j(qi)B2cr{qi)v2 . 

Dependencies  of  critical  magnetic  inductions  upon  q/n  obtained  in  different  approaches  for 
various  values  of  flight  Mach  number  are  shown  in  the  Fig.7.  MHD  generator  is  located  in 
cross-section  2  of  scramjet.  It  is  easy  to  see  that  critical  magnetic  induction  is  monotone 
increasing  function  of  ratio  q/n2 .  Increase  of  flight  Mach  number  leads  to  decrease  of  the 
critical  magnetic  induction.  The  critical  magnetic  inductions  calculated  both  in  the  second 
and  in  the  third  approaches  are  less  than  one  for  the  first  approach.  In  the  range  of 
0.\<q/n<\Q  the  deviations  for  different  approaches  are  insignificant.  In  the  range  of  small 
values  of  power  spent  on  flow  ionization  (q/n2< 0.01)  the  second  approach  belittles  value  of 
critical  magnetic  induction  noticeably.  It  is  consequence  of  electron  attachment  processes, 
which  cannot  be  taken  into  account  in  the  second  approach  in  principle.  Nevertheless  the 
second  approach  provides  good  agreement  with  numerical  results  obtained  in  the  first 
approach  in  wide  range  of  parameter  q/n2  variation.  Besides,  in  the  second  approach  we 
obtain  simple  analytical  formula  for  critical  value  of  magnetic  induction  : 

Bcr  =  ^menkc  Jq~W~P /(/c  (l  -  /c  )  eV)  (12a) 


This  formula  gives  us  opportunity  to  understand  dependence  of  the  critical  magnetic 
induction  on  the  pivotal  parameters.  At  fixed  value  of  magnetic  induction  we  will  use  another 
notion  namely  the  critical  power  density  spent  on  flow  ionization  with  corresponding  notation 
qcr.  In  the  second  approach  it  can  be  obtained  by  converting  expression  (12a). 


1  rk(l-k)e2ff2v2Y 
Wip  \  mekcn  , 


(12b) 


Self-sustained  operational  mode  of  MHD  generator  with  non-equilibrium  conductivity  is 
realized  then  0 <qi<qcr- 

Developed  model  allows  us  to  formulate  requirements  for  parameters  of  ionizer  and 
MHD  generator  at  which  MHD  generator  with  nonequilibrium  conductivity  will  operate  in 
self-sustained  mode.  Moreover  we  can  suppose  algorithms  for  optimal  choosing  of  the  MHD 
parameters  by  using  the  discussed  results.  The  next  simple  analysis  demonstrates  the 
possibilities  for  optimal  choice  of  ionizer  and  magnetic  system  parameters. 

MHD  generator  with  non-equilibrium  conductivity  is  a  part  of  complex  gas-dynamic 
system,  so  the  electric  energy  production  is  not  the  only  its  function.  MHD  generator,  located 
in  scramjet,  must  effectively  regulate  characteristics  of  scramjet.  According  to  next  chapter, 
to  noticeably  improve  scramjet  performance  it  is  necessary  to  ensure  MHD  interaction 
parameter  Sv>0.1,  where  Sv  =  <jB2  Lf  pv ,  L  is  the  length  of  MHD  generator.  Decrease  the 


power  density  input  to  ionization  causes  the  flow  conductivity  to  be  decreased.  Thus  to 
ensure  required  MHD  interaction  parameter  while  decreasing  qi  it  is  necessary  to  increase 
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integral  quantities:  Wion  is  the  power  put  to  flow  ionization  in  MHD  generator  channel,  Wg  is 
the  power  produced  by  the  MHD  generator.  The  power  Wg  is  determined  by  the  ratio 

L 

w,=-\o  •E)A(x)dx ,  where  L  is  the  MHD  generator  length.  The  power  Wion  is  determined 

o 

L 

by  the  ratio  Wion  =  J  qiA(x)dx .  Also  we  introduce  factor  r  as  the  ratio  of  power  input  to 

0 

ionization  and  power  produced  by  MHD  generator  r=Wion/JVg.  It  is  obvious  that  self- 
sustained  operational  mode  is  realized  when  0<r<l. 

The  power  Wg  can  be  easily  calculated  in  terms  of  the  enthalpy  extraction  ratio  r\: 
Wg=rj-W0,  where  W0  =  Ax px y cp Txa>(Mx) .  To  calculate  Wion  we  take  into  account  that 

A(x)=AipiVi/p(x)v(x).  Change  of  variables  in  integrating  over  x  by  using  (17,19)  let  us 
obtain  next  formula  for  factor  r: 


i  n 

iv)=—\- 


drj 


77  J0  k(\-k)x<jB2v2  /qt  ^  ^ 

Let  us  assume  that  quantity  qion  is  constant  in  the  whole  of  MHD  generator  volume.  By 
substituting  (17,  20)  in  (22)  we  obtain: 

f  .  v-6 


*(//)  = 


\dcr) 


1  '' 

-I 


t\2bD 


Z(k,  y/,Ml,r/') 


-drj' 


In  the  case  with  a  constant  e-beam  current  density  we  obtain: 


•(/;)  = 


d 

WcrJ 


(1 +b)D 


1 


-dr/' 


Here  the  local  critical  power  density  qcr  is  determined  by  the  ratio: 


da  h  =  dobk (1  -  k)  JO-0  (p0  /  px  f  B 


>2.  .2 


(23a) 


(23b) 


(24) 


When  we  use  the  second  approach  in  determining  of  nonequilibrium  conductivity  the  local 
critical  power  will  be  calculated  under  the  formula  (12b).  Formulae  (23)  determine 
requirements  to  ionizer  and  MHD  generator  to  maintain  self-sustained  operational  mode. 

Fig  9  shows  function  r( q)  for  Faraday  and  Hall  MHD  generators  at  various  values  of 
Mach  number  Mi  and  various  values  of  a  load  factor  k.  In  all  the  cases  MHD  flow  is 
characterized  by  the  value  £  =  |(^) .  One  can  see  that  function  r(rj)  is  monotone  increasing 

function.  Value  7=77*  defined  from  relation  r(rf)=\  determines  upper  limit  of  the  enthalpy 
extraction  ratio  at  which  self-sustained  operational  mode  of  MHD  generator  with 
nonequilibrium  conductivity  exists.  One  can  see  that  values  7*  for  Faraday  MHD  generator  is 
greater  than  for  Hall  MHD  generator  at  the  same  conditions.  It  follows  from  Fig.  9  that  range 
of  existence  of  the  self-sustained  operational  mode  is  extending  both  for  Faraday  and  Hall 
MHD  generators  while  increasing  Mach  number.  Increase  of  load  factor  leads  to  extension  of 
the  range  of  self-sustained  operational  mode  for  Faraday  MHD  generator.  For  Hall  MHD 
generator  the  same  effect  is  achieved  by  decreasing  the  load  factor. 

Fig.  10  demonstrates  relative  power  characteristics  both  for  Faraday  and  Hall  MHD 
generators.  Here  and  hereinafter  to  clarity  we  use  the  designation  qg  instead  of  7  to  underline 
that  7 g=Wg/Wo  is  relative  power  produced  by  MHD  generator.  The  factor  7 ion=Wion/Wo  is  the 
relative  power  spent  on  flow  ionization.  The  difference  r/g-r}ion=(Wg-Wion)/Wo  determines  the 
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Fig.2.  Electron  beam  power  loss  in  various  cross-sections  of  scramjet  as  a  function  of  e-beam 
current  density,  £6=100  keV.  Numbers,  which  denote  the  curves,  correspond  to  location  of 
the  cross-sections  according  to  Fig.  1 .  Flight  Mach  number  Mx=  8,  free-stream  dynamic 
pressure  is  equal  to  40  kPa. 
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pvA  =  m  =  const 

dv  dp 
pv —  +  — =  0 
dx  dx 


f 


pv 


dx 


.2  \ 


CPT+— 

KP  2y 


=  q(x)  +  qr 


(27) 


p  =  RpT 

where  q(x)  is  power  density  released  in  combustion  chamber  in  result  of  fuel  combustion,  qr 
is  density  of  heat  release  in  result  of  plasma  recombination  processes.  Total  heat  addition  in 

L  L 

combustion  chamber  Q=Qj+Qr,  where  Qf  =  J  q  (x)  A  (x)  dx ,  Qr  =  jqr- A  (x)  dx  =  Wion . 

o  o 

Let's  consider  the  combustion  chamber  working  in  a  mode  with  a  constant  pressure 
and  in  a  mode  with  a  constant  density.  It  is  evidently  from  the  set  (27)  that  flow  parameters 
modification  in  combustion  chamber  with  a  constant  pressure  are  determined  by  the  ratios: 


r3  =t2+at+ 


w:. 


mc_ 


P3=Pl’  V3  =V2 

^  H: 
AT  = 


(28a) 


cp(aL0+l) 

where  Hu  -  a  calorific  value  of  fuel,  Lo  -  stochiometric  factor,  a  is  air-fuel  ratio  (a>l). 

Corresponding  ratios  for  combustion  chamber  with  a  constant  density  have  a  form: 

/  \ 

W;„ 


Ts-Ti+y' 


AT  +  - 


me 


p  J 


P3  =  Pi  ■  (T3 1  Tl) 
vl=vl.(i-i/ryCpiT3-T2) 


(28b) 


3-4.  MHD  accelerator 

The  MHD  accelerator  is  characterized  by  parameters  £,3,  kj.  (For  MHD  accelerator  the 
load  factor  ^3  is  greater  than  unity  ks>l)  It  is  supposed  that  part  (1-r)  of  power  produced  by 
MHD  generator  is  transferred  to  MHD  accelerator  Wa  =  (l-r)lfg  .  Flow  parameters  in  MHD 

accelerator  are  described  according  to  theory  developed  in  [8,  9].  The  change  of  flow 
parameters  in  a  channel  of  MHD  accelerator  is  determined  as  follows: 

7;=r3+^(l  +  ^)T„-®(M0)-I7-(l-r) 

/V3 


1L  =  2L  + 


v;  1  -  £3  •  (k3  - 1) 


cT0-G>(M0)-Tp(l-r) 


tv  2 

£± 

P3 


fr\D^ 


\Tu 


(29) 
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where  A  = 


r_ _ 4 

Y- 1^+1 


—  1. 


4-5.  Nozzle 


We  suppose  that  nozzle  flow  is  isentropic.  Thus  the  relationship  between  a  relative 
change  of  a  flow  pressure  and  a  relative  change  of  a  flow  temperature  in  the  nozzle  is  as 
follows: 


Ps_  = 
Pa 


2 


\Ta  J 


r- 1 


y+c,-(r.-rs) 


(30) 


We  consider  nozzle  operating  at  design  conditions,  so  the  set  of  equations  (25,26,28- 
30)  can  be  closed  in  supposing  the  pressure  at  an  exit  of  the  nozzle  is  equal  to  pressure  in  a 
surrounding  medium  p5=po ■  In  this  situation  we  use  evident  relation: 


Pl  =  Pl.p±.Pl.H.Pl  =  i 
Po  Pa  A  Pi  A  Po 


In  view  of  this  ratio,  we  can  obtain  formula  for  calculation  of  a  flow  temperature  at 
the  exit  of  the  nozzle  from  the  set  of  equations  (25,26,28-30).  In  the  case  when  combustion 
chamber  working  in  a  mode  with  a  constant  pressure  we  obtain: 


Ts  = 


(i-i  ir) 


Tx 


U  J 


+i 


f  rr,  \ 


V^3  , 


#3 

#3+1 


(31a) 


In  the  case  when  combustion  chamber  working  in  a  mode  with  a  constant  density  we 

obtain: 


t5  = 


Vi+1  f  T,  V1  Vr)(T.  V3+1 


(31b) 


V"i  J 


l2  J 


V^3  J 


The  exhaust  velocity  of  gas  from  the  nozzle  is  determined  through  temperature  A  with 
use  of  an  energy  conservation  law  by  the  following  ratio: 

vs = A2+2cA+Ar-ri)  <32> 


The  obtained  ratios  allow  calculating  specific  impulse  of  scramjet  with  MHD  control 
Isp.  An  explicit  formula  for  the  specific  impulse  of  air-breathing  propulsion  systems  when  a 
mass  flow  rate  of  a  fuel  is  neglected  in  a  comparison  with  a  mass  flow  rate  of  an  air, 
according  to  [1 1]  is  as  follows: 

rvT 

ISP=—  (AA5-v0)  (33) 

g 


where  (Pn  is  the  nozzle  non-ideality  factor  (in  calculations  we  suppose  that  q>N=  I ). 

The  set  of  the  formulas  (25,26,28-33)  allows  one  calculating  specific  impulse  of  the 
scramjet  with  MHD  control  at  given  parameters  of  the  inlet,  MHD  systems  and  the 


26 


combustion  chamber.  To  determine  the  range  of  parameters  at  which  MHD  interaction 
improves  scramjet  performance  we  use  an  obvious  functional  ratio  (d  Isp/<r?r/^  q  >  0 . 

Having  done  necessary  transformations  we  obtain  requirements  for  £/  parameter.  In 
the  case  of  combustion  chamber  working  in  a  mode  with  a  constant  pressure,  MHD 
interaction  increases  specific  impulse  of  scramjet  if  the  next  inequality  is  true: 

tf-(‘T(1-r(0)))  <34a) 

It  is  evidently  that  r(0)  can  be  determined  as  the  ratio  of  power  densities  r(0)=qi  /qg.  As 
0<r(0)<l  and  kj>l  the  factor  8  has  a  value  in  the  range  0<8<1.  Inequality  (34a)  combines 
some  parameters  of  scramjet,  MHD  systems  and  ionizer.  Cursory  analysis  of  the  inequality 
let  us  to  do  the  inference  that  internal  MHD  generator  allows  one  to  increase  the  scramjet 
specific  impulse  only  when  we  use  MHD  generator  in  regime  with  pressure  increasing  along 
the  channel  length.  Indeed  for  MHD  generator  y/j>0  and  0<8<1,  so  the  right  part  of  (34a)  is 
positive.  According  to  [8],  condition  %i>0  corresponds  to  MHD  flow  with  dp/dx>0.  If  we 
consider  MHD  generator  with  constant  cross-sectional  area  the  inequality  (34a)  takes  on  the 
next  form: 


y  +  l/V!  Tx 
y  —  1  +  - — 

M,2- 1  AT 


i+A 


(34b) 


This  inequality  prescribes  the  range  of  variation  of  MPCE  subsystems  parameters  at  which 
using  of  MHD  generator  with  a  constant  cross-sectional  area  allows  one  to  increase  scramjet 
specific  impulse. 

According  to  Fig.  16  and  17  the  temperature  77  depends  on  the  flight  Mach  number 
Mo.  The  value  Mi  in  (34b)  according  to  (26)  depends  on  the  flight  Mach  number  too.  Thus  in 
result  of  analysis  of  the  inequality  (34b)  we  obtain  limitations  on  flight  Mach  number  at 
which  MHD  interaction  causes  the  MPCE  specific  impulse  to  be  increased.  Fig.  18  shows 
dependencies  of  maximal  flight  Mach  number  upon  factor  5  for  various  values  of  y/  factor. 
MHD  interaction  increase  specific  impulse  of  scramjet  in  the  case  when  flight  Mach  number 
Mao  <Mmax.  As  it  follows  from  (34a)  factor <7 mainly  characterize  relative  power  spent  on  flow 
ionization.  In  situation  when  power  spent  on  ionization  is  equal  to  zero  (initially  ionized 
flow)  8=0.  Factor  8=1  describes  situation  for  which  whole  of  power  produced  by  MHD 
generator  is  spent  on  flow  ionization.  It  is  evidently  that  while  increasing  the  factor  8  the 
Mmax  is  decreasing.  Decrease  of  the  total  turning  angle  On  causes  the  maximal  Mach  number 
Mmax  to  be  increased.  At  small  value  of  8,  minimization  of  y/  factor  provides  expansion  of 
range  of  Mach  number  at  which  MHD  interaction  improves  scramjet  performance.  Contrary 
at  S~1  increase  of  ^factor  causes  the  Mmax  to  be  increased. 

In  the  case  of  combustion  chamber  working  in  a  mode  with  a  constant  density,  MHD 
interaction  increases  specific  impulse  of  scramjet  if  the  next  inequality  is  true: 


6>(r-i) 


AT 


(35a) 


If  we  consider  MHD  generator  with  constant  cross-sectional  area  the  inequality  (35a)  takes 
on  the  next  form: 
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MHD  interaction  parameters.  For  Faraday  MHD  generator  increase  of  q/qCr  ratio  leads  to 
decrease  of  optimal  values  both  for  ki  factor  and  MHD  interaction  parameter.  So  we  can 
conclude  that  ratio  q/qcr  is  important  parameter  that  influence  on  optimal  value  both  for  load 
factor  and  MHD  interaction  parameter.  According  to  previous  chapter  this  ratio  is  determined 
by  parameters  of  ionizer,  magnetic  system  and  flow  parameters.  To  formulate  requirements  to 
these  subsystems  we  have  calculated  characteristics  of  scramjet  using  dimensional 
parameters. 

Fig.31  shows  contours  for  relative  specific  impulse  and  power  excess  produced  by 
MHD  generator  in  scramjet  with  MHD  control  in  the  axes:  magnetic  induction  and  power 
density  spent  on  ionization.  These  results  are  obtained  for  free  stream  dynamic  pressure  40 
kPa.  These  dependencies  correlate  with  the  qualitative  results  presented  in  the  Fig. 8.  Namely 
there  are  optimal  values  of  power  density  qf  at  which  required  both  specific  impulse  and 
power  excess  can  be  achieved  at  minimal  value  of  magnetic  induction. 

As  it  follows  from  the  Fig.32-33  optimal  value  of  the  load  factor  for  specific  impulse 
differ  from  one  for  power  excess.  Maximal  power  excess  is  achieved  at  ki  near  0.5  both  for 
case  with  given  power  density  spent  on  ionization  (Fig.32)  and  case  with  given  magnetic 
induction  (Fig.33).  As  for  specific  impulse  its  maximal  values  can  be  achieved  at  the  load 
factors,  which  are  essentially  less  than  0.5. 

Obtained  results  show  that  MHD  generator  with  nonequilibrium  conductivity  can  be 
used  on  hypersonic  aircraft  both  for  electric  power  production  and  for  scramjet  specific 
impulse  increasing.  Requirements  for  parameters  of  the  MHD  generator  depend  on  its 
functionality. 


Fig.  13  Simplified  scheme  of  scramjet  with  MHD  control  under  “AJAX”  concept,  according 
to  [11] 
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Fig.  14  Simplified  scheme  of  scramjet  with  MHD  control  in  quasi-l-D  model 

0-1  -  inlet;  1-2  MHD  generator;  2-3  -  combustion  chamber;  3-4  -  MHD  accelerator;  4-5  - 

nozzle 


Fig.15.  Density  contours  in  scramjet  inlet.  Md=10,  0n=15°,  Fth=0.12,  Mo=10 
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Fig.20  Relative  specific  impulse  (a)  and  relative  power  spent  on  ionization  (b)  for  MHD 
generator  with  various  parameters  £/•  6n=\ 5°,  Mj=\ 0,  F,/,=0. 1 2,  M0= 6,  q/qcr=0.05,  power 
density  is  constant  in  MHD  generator  volume. 
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Fig.22  Relative  specific  impulse  vs  load  factor  for  scramjet  with  Faraday  MHD  generator  for 
various  values  of  enthalpy  extraction  ratio.  #y=15°,  Mj=  10,  Fth=0. 1 7,  Mo=6,  qi/qcr= 0.01, 
k3= 1.1 
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b) 

Fig.24  Relative  specific  impulse  (a)  and  power  produced  by  MHD  generator  in  scramjet  vs 
MHD  interaction  parameter  for  various  relative  values  of  inlet  throat  Fth .  $v=15°,  M(/=10, 
Mq=6,  q/qa-=0,  kf=  0.5,  kj=\.2. 
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4.  TWO-DIMENSIONAL  MODEL  FOR  MHD  CONTROLLED  INLET  OF  SCRAMJET 


This  chapter  is  devoted  to  development  of  mathematical  model  of  combined  gas- 
dynamical  system  which  includes  MHD  generator  and  ionizer  in  two-dimensional  approach 
and  numerical  calculations  of  flow-fields  in  gas-dynamical  channel  in  a  wide  spectrum  of 
flow  parameters,  configuration  of  gas-dynamical  system,  power  spent  on  flow  ionization, 
shape  and  length  of  plasma  formation.  As  the  combined  gas-dynamical  system  we  consider 
scramjet,  which  includes  MHD  generator  with  non-equilibrium  conductivity,  located  in  inlet 
of  the  scramjet.  It  is  supposed  that  MHD  generator  allows  us  to  control  flow-field  and 
improve  scramjet  performance. 

Typical  geometry  of  scramjet  inlet  with  design  Mach  number  Mj=  10,  total  turning 
angle  $v=15°  and  relative  value  of  inlet  throat  Fth= 0.1  is  shown  in  the  Fig.34  ( Fth=l-Yb/Ya ). 
Arbitrary  ionized  region,  where  MHD  interaction  is  effected,  is  shown  in  the  Fig.34.  We 
consider  homogeneous  magnetic  field  which  has  both  x  and  y  components  (B  =  {BX,B  , o}). 


Magnetic  field  will  be  characterized  by  two  parameters:  modulus  of  magnetic  induction 
=M  +  B]  and  angle  of  magnetic  field  orientation  a  ( Bx  =  B  ■  cos(a) , 


B  = 


B 


By  =  B  ■  sin  (a) ).  It  is  supposed  that  electrodes  of  MHD  generator  are  located  in  the  plane  of 

the  figure  at  coordinates  z=-co  and  z=oo.  Conductivity  of  flow  is  created  by  using  of  e-beam 
ionizer.  Flow  conductivity  depends  upon  the  power  density  spent  on  flow  ionization  and  flow 
parameters.  To  calculate  flow  conductivity  we  use  formulas  presented  in  the  second  chapter. 
Position  and  shape  of  ionized  region  are  chosen  arbitrarily  to  investigate  how  the  geometrical 
characteristics  of  the  ionized  region  influence  on  the  inlet  performance. 

Flow- field  in  the  MHD-controlled  inlet  of  scramjet  is  calculated  in  two-dimensional 
Euler  approach.  Set  of  equations  used  in  calculations  is  listed  below: 

dx  dy 


d(pv1  +  p)  |  dpvxvy 

dx  dy  x’  (36) 

d(pv:+p)  dpvxvy  _ 
dy  dx  Jy ’ 


d{e  +  p)vx  ^  d{e  +  p)vy 
dx  dy 

e  =  P(cvT  +  (v;  +  vl)l2\ 


-  Qg + qr  > 
p  =  RpT , 


where  p  is  the  flow  density,  vx  and  vy  are  the  x  and  y  components  of  flow  velocity,  p  is  the 
static  pressure,  T  is  the  flow  temperature,  cv  is  the  specific  heat,  R  is  the  absolute  gas 
constant.  Expressions  for  the  right  side  of  equations  (36)  have  the  following 
view:  f  =  {fx,fy,  o},  f  =  jxB,  qg=j-E,  where  j  is  the  current  density,  E  is  the  electric  field. 

The  generalized  Ohm’s  law  determines  relations  between  electromagnetic  components: 

y+//(yxs)  =  (7(f?+vxi?)  (37) 


where  p  is  the  electron  mobility,  cr  is  the  electrical  conductivity.  The  quantity  q,  included  in  a 
right  side  of  (36)  determines  the  power  density  put  into  a  flow  as  a  result  of  a  recombination 
of  charged  particles. 


48 


temperature  behind  the  shock  when  oblique  shock  is  deflected  by  magnetic  field  it  is 
necessary  to  use  magnetic  field  with  a=90°.  Otherwise  it  is  necessary  to  apply  magnetic  field 
with  a<90°. 

Now  we  consider  peculiarities  of  MHD  interaction  on  stream  flowing  around  the 
wedge  when  a>90°.  Pressure  colour  map  in  the  flow  around  the  wedge  with  MHD  control 
with  magnetic  field  which  is  characterized  by  a=135°  is  shown  in  the  Fig.42.  In  this  case  the 
Lorentz  force  is  negative  and  we  expect  that  MHD  interaction  decrease  the  angle  of  oblique 
shock.  In  comparing  the  Fig.42  and  Fig.39  one  can  conclude  that  oblique  shock  position  in 
this  case  is  practically  similar  to  position  of  oblique  shock  without  MHD  interaction.  It  can  be 
explained  by  pressure  increase  behind  the  shock.  Fig.43  demonstrates  that  pressure  increase 
behind  the  shock  in  this  case  is  significant  in  comparison  with  a=45°  MHD  interaction 
(Fig.40a)  and  practically  the  same  as  in  the  case  of  MHD  interaction  with  a=90°  (Fig.40b). 
According  to  the  Fig.44  the  temperature  increase  behind  the  shock  for  a=135°  is  essentially 
less  when  one  for  a=90°  (Fig.41b).  Fig.45  demonstrates  positions  of  oblique  shocks  in  the 
case  of  MHD  interaction  with  a=45°  and  a=135°  and  in  the  absence  of  MHD  interaction. 
One  can  see  that  MHD  interaction  with  a=45°  increases  angle  of  oblique  shock  and  MHD 
interaction  with  a=135°  very  slightly  decreases  angle  of  oblique  shock.  Such  effect  is 
obtained  for  5=5T  (a=135°)  and  B=  IT  (a=45°).  So  one  can  conclude  that  decreasing  the 
angle  of  oblique  shock  by  MHD  control  is  more  complicated  problem  than  increasing  one. 

And  now  we  consider  possibilities  of  MHD  control  for  modification  of  flow  field  in 
inlet  at  off-design  conditions.  Firstly  situation  for  flight  Mach  number  greater  than  design 
Mach  number  will  be  investigated.  In  papers  [12,13]  it  was  shown  that  MHD  control  in  this 
case  allows  us  to  modify  flow  field  in  such  a  way  that  it  becomes  like  a  flow  field  in  design 
conditions.  Thus  in  this  report  we  only  briefly  consider  how  various  plasma  configurations 
influence  on  the  flow  field  at  Mo>Md. 

The  Fig.46  demonstrates  possibilities  of  MHD  interaction,  which  is  located  in  the 
range  0<x<0.5,  to  modify  flow  field  in  scramjet  inlet.  Fig.46a  shows  Mach  number  contours 
in  inlet  at  off-design  conditions  without  MHD  interaction  at  flight  Mach  number  Mo=12.  In 
this  case,  unlike  to  design  conditions  (see  the  Fig.35),  the  oblique  shocks  intersection  occurs 
not  in  cowl  lip.  Increase  of  power  density  spent  on  ionization  (Fig.46  b-d)  causes  the  oblique 
shocks  angles  to  be  increased.  When  qi=0.022W/cms  oblique  shocks  are  intersected  in  the 
cowl  lip  and  flow  field  becomes  quite  similar  to  the  one  shown  in  the  Fig.35.  When  qt 
becomes  greater  than  0.022  W/cm3  (Fig.46d)  the  flow  field  becomes  similar  to  situation  with 
the  flight  Mach  number  is  less  than  design  one. 

In  the  Fig.47a  MHD  interaction  is  realized  in  the  range  3<x<3.5.  In  this  situation 
oblique  shocks  intersection  occurs  not  exactly  in  cowl  lip.  But  the  point  of  shocks 
intersection  is  nearer  to  cowl  lip  than  in  the  case  without  MHD  interaction  (Fig.46a).  In  the 
Figs.47b-c  MHD  interaction  is  realized  in  the  range  with  3<x<4  and  y<0.7,  which  is 
schematically  shown  in  the  Fig  37b.  When  magnetic  field  is  orthogonal  to  x-axis  (a=90°)  the 
shocks  intersection  occurs  not  exactly  in  cowl  lip  (Fig.47b).  But  in  changing  the  magnetic 
field  orientation  (a=45°)  point  of  oblique  shocks  intersection  moves  to  cowl  lip  (Fig.47c). 

Fig.47d  corresponds  to  ionized  range  located  in  0<x<5  and  y<0.2-x,  which  is 
schematically  shown  in  the  Fig  37d.  One  can  see  that  in  this  case  the  flow  field  is  quite 
similar  to  the  flow  field  at  design  conditions.  Fig.48  demonstrates  density  distributions  in 
cross-section  of  MHD  controlled  inlet  at  x=7.1  for  the  configuration  of  ionized  region.  It 
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normalized  flow  compression  (tt/tto),  where  n =(pi/p0)  is  the  flow  compression  in  MHD 
controlled  inlet  and  Tio  is  the  flow  compression  in  the  inlet  without  MHD  control.  Fig.  54 
demonstrates  how  the  normalized  relative  ratios  depend  on  the  orientation  angle  of  the 
magnetic  induction  a  for  various  values  of  B.  Flight  Mach  number  Mo=8,  Md=10.  Ionized 
region  is  a  band  located  in  3<>c<4.  One  can  see  that  both  for  mass  flow  and  for  flow 
compression  there  are  optimal  values  of  a  angle,  greater  than  90°,  at  which  these  parameters 
achieve  a  maximum.  Maximal  air  capture  in  inlet  increase  while  increasing  the  magnetic 
induction  from  1  to  3  Tesla. 

Fig.  5  5  shows  dependencies  of  cp/cpo  as  a  function  of  the  load  factor  k  for  various 
orientations  of  magnetic  field  at  value  B=3T.  Results  of  calculations  correspond  to  flight 
Mach  number  Mo= 8.  Ionized  region  is  a  band  located  in  3<x,<4.  It  is  easy  to  see  that  in  the 
considered  situation  MHD  control  allows  us  to  increase  the  air  capture.  There  are  some 
optimal  values  for  load  factor  (depending  on  the  angle  a)  at  which  maximal  air  mass  flow  in 
MHD  controlled  inlet  can  be  achieved.  Maximal  air  capture  in  inlet  increase  while  increasing 
the  a  angle  from  130°  to  145°.  Fig.56  demonstrates  dependencies  of  normalized  flow 
compression  n/no  as  a  function  of  the  load  factor  k  in  conditions  similar  to  the  Fig.56.  One 
can  see  that  flow  compression  in  the  inlet  monotonically  decrease  while  increasing  the  load 
factor.  Load  factor  optimal  for  increasing  the  air  mass  flow  is  not  optimal  for  flow 
compression  increasing.  Fig.57  and  Fig.58  demonstrate  dependencies  of  cp/cpo  and  tc/tco, 
correspondingly,  as  a  functions  of  the  load  factor  k  for  various  a  angle  at  value  B=3T  for 
flight  Mach  number  Mo=6.  These  dependencies  are  similar  to  corresponding  dependencies 
obtained  for  Mo= 8  (Fig.55,56).  Thus  obtained  results  show  that  at  given  region  of  MHD 
interaction  in  inlet  the  air  capture  and  the  flow  compression  noticeably  depend  on  the  value 
and  orientation  of  magnetic  induction  and  the  value  of  the  load  factor. 

Now  we  consider  influence  of  shape  and  position  of  ionized  region  on  characteristics 
of  the  MHD  controlled  inlet.  Fig.59  shows  dependencies  of  cp/cpo  as  a  function  of  position  of 
ionized  region  Xi  in  x-direction  for  various  height  of  the  region  17.  Ionized  region  is  column, 
shown  in  the  Fig.37b,  located  in  X/<x<Xi  +  l  and  yXY2.  Air  mass  flow  monotonically 
increases  while  increasing  the  height  of  ionized  region.  For  Y2>1.25  the  air  mass  flow 
practically  doesn’t  depend  on  the  height  and  coincides  this  the  one  for  ionized  region  in  form 
of  band.  Dependencies  of  cp/cpo  upon  position  Xt  is  more  complicated.  At  small  heights  Y2  the 
air  mass  flow  monotonically  increases  while  increasing  the  Xj  value.  For  Y2>1.1  there  are 
optimal  positions  of  Xi  at  which  maximal  air  mass  flow  is  achieved.  The  optimal  value  for  Xi 
decreases  while  increasing  the  height  of  the  ionized  region.  Fig.60  shows  dependencies  of 
n/no  as  a  function  of  position  of  ionized  region  Xj  for  various  heights  of  the  region  Y2. 
Conditions  are  the  same  as  in  the  Fig.59.  Similarly  to  the  air  mass  flow  in  inlet  the  flow 
compression  increases  while  increasing  the  height  of  the  ionized  region.  Positions  of  ionized 
region  Xj  optimal  for  increasing  the  flow  compression  are  slightly  differ  from  positions 
optimal  for  increasing  the  air  mass  flow. 

Results  of  investigation  of  MHD  controlled  inlet  in  the  case  when  ionized  region  has  a 
rectangular  shape:  3<x<4  and  Yt<y<Y2  (shown  in  the  Fig.37c)  are  presented  in  the  Fig.61  and 
62.  The  Fig.61  shows  ratio  cp/cpo  as  a  function  of  bottom  border  coordinate  Yj  at  various 
coordinates  of  the  upper  border  Y2  for  the  ionized  region.  One  can  see  that  the  air  mass  flow 
in  the  inlet  monotonically  increases  while  increasing  the  Y2  value.  For  each  position  of  Y2 
there  is  optimal  value  of  Yi,  which  is  greater  than  0.8m  in  considered  configuration.  Fig.62 
shows  dependencies  of  n/no  as  a  function  of  position  of  ionized  region  Yt  for  various 
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Fig.37  Simple  shapes  of  ionized  region  used  in  calculations  of  flow-fields  in  MHD-controlled 
inlet 


X 


Fig.42  Pressure  colors  map  in  stream  flowing  the  wedge  with  0=10°.  Mo=8,  t/,=4W/cm3, 
A=0.5,  B=3T,  a=135°.  Ionized  region  located  in  2<x<2.5. 


Fig.43  Pressure  contours  in  stream  flowing  the  wedge  with  0=10°.  Mo=8,  g;=4W/cm  ,  k= 0.5. 
B=3T,  a=135°.  Ionized  region  located  in  2<x<2.5. 
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b)  q, .  w/cm3 

Fig.50.  Relative  mass  flow  rate  and  pressure  increase  in  MHD  controlled  inlet  vs  power 
density  spent  on  ionization  for  various  magnetic  field  orientations.  Ionized  region  is  located 
in  xi<x <ki+1.Mo=8,  Mj=  10,  6n=  15°,  F,/;=0.12,  5=1  T,  k= 0.5,  a  values  are  shown  in  the 
figure,  (a)  x/=0,  (b)  x/=l . 
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Fig.56.  Normalized  relative  pressure  increase  in  MHD  controlled  inlet  vs  load  factor  for 
various  a.  Mo= 8,  Md=  10,  6!v=15°,  F,/,=0.12,  B= 3T,  </,=l  W/crn3.  Ionized  region  is  located  in 
3<x<4. 


Fig.57.  Normalized  relative  air  mass  flow  in  MHD  controlled  inlet  vs  load  factor  for  various 
a.  Mq= 6,  Mj=  10,  6n=  15°,  F,/,=0. 1 2,  B= 3T,  qr\  W/cm3.  Ionized  region  is  located  in  3<x<4. 
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conductivity  of  flow.  Also  we  will  consider  simple  model  of  low  ionized  plasma  consisting  of 
neutral  molecules,  electrons,  negative  and  positive  ions  with  corresponding  number  densities: 
n,  ne,  n.,  n+.  Set  of  kinetic  equations  for  the  plasma  components  concentration  is  the  next: 

Qyi 

~^-  =  I  +  kd-n-ne- kaN20ne  -  (3n+ne 

^-  =  1-  pun+n_-  pn+ne  (1) 

ct 
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=  kaNl2ne  ~  KlN02  ■  «-  ~  Pu>hn 


where  /  is  e-beam  induced  ionization  rate,  N0i  -  is  concentration  of  oxygen  molecules  in  air, 


ka  and  kd  are  the  rate  constants  of  attachment  and  detachment  of  electrons  respectively,  [5  and 
Pa  are  the  rate  constants  for  electron-ion  and  ion-ion  recombination  respectively.  These 
constants  are  functions  of  gas  temperature  Tg  and  electron  temperature  Te.  In  calculations  we 
use  dependencies  from  papers  [4,  5],  Basic  reactions  and  corresponding  rate  coefficients  are 
shown  below: 
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Ionization  rate  /  is  determined  in  terms  of  e-beam  characteristics  by  the  ratio: 

i  =  3l  Ah!e)pY(Eb) 

W,  W; 


(6) 


where  jb  -  is  the  e-beam  current  density,  Eb  -  is  the  energy  of  electron,  Y  is  the  electron 
stopping  power, 


magnetic  induction  value.  The  decreasing  curve  in  Fig.  8  corresponds  to  constant  value  of 
MHD  interaction  parameter  5V=0.1.  The  increasing  curve  in  Fig.8  corresponds  to  critical 
regime  of  MHD  generator  for  which  qg=qi.  Parameters  B  and  qt  which  are  located  in  shaded 
area  correspond  to  self-sustained  operational  mode  of  MHD  generator  with  non-equilibrium 
conductivity  with  MHD  interaction  parameter  5V>0. 1 .  It  follows  from  Fig.8  that  there  is 
optimal  value  of  energy  input  to  ionization,  which  ensures  required  regime  of  MHD  generator 
in  scramjet  at  minimal  magnitude  of  magnetic  induction. 

Table  2 

Optimal  values  of  the  power  spent  on  flow  ionization  qopt  and  the  magnetic  induction  Bopt  for 
scramjet  with  various  values  of  the  total  turning  angle  On  for  various  values  of  flight  Mach 
number  Mo  at  free  stream  dynamic  pressure  40kPa,  S*=0.1,  L=lm. 


On, 

degree 

Mo 

Position  of  MHD  generator  in  scramjet  according  to  the  Fig.l 

] 

2 

3 

Qoph 

W/cm 3 

Bmin,  T 

tfopti 

W/cm3 

Bmin,  T 

Qopb 

W/cm 3 

Bmm,  T 

10 

6 

5.66 

1.1 

8.61 

1.5 

16.2 

2.7 

8 

9.04 

0.75 

15.3 

1.1 

29.7 

2.0 

10 

12.1 

0.55 

24.3 

0.94 

47.5 

1.7 

12 

16.4 

0.45 

35.2 

0.80 

58.0 

1.2 

15 

6 

6.35 

1.2 

11.2 

1.9 

14.3 

2.7 

8 

9.94 

0.83 

20.5 

1.5 

26.8 

2.1 

10 

14.3 

0.63 

34.9 

1.3 

44.4 

1.7 

12 

21.0 

0.54 

53.8 

1.1 

54.4 

1.3 

20 

6 

7.35 

1.3 

14.4 

2.4 

12.0 

2.5 

8 

12.0 

0.96 

27.7 

2.0 

23.0 

2.2 

10 

17.0 

0.73 

45.5 

1.6 

37.9 

1.8 

12 

23.1 

0.59 

67.2 

1.4 

44.5 

1.3 

Optimal  values  of  power  density  spent  on  flow  ionization  and  corresponding 
minimum  magnitudes  of  magnetic  induction  are  presented  in  the  Table2  for  various  flight 
Mach  numbers.  Calculations  were  made  for  different  inlets  with  the  total  turning  angle 
varying  from  10  to  20  degree.  One  can  see  that  minimal  magnetic  induction  decreases  while 
increasing  the  flight  Mach  number.  Magnetic  induction  in  internal  part  of  inlet  needs  to  be 
greater  than  in  external  one.  The  total  turning  angle  increase  causes  both  the  minimal 
magnetic  induction  and  the  optimal  power  density  spent  on  flow  ionization  (in  external  part 
of  inlet)  to  be  increased.  The  qopt  value  in  internal  part  of  inlet  is  greater  than  in  external  one. 
The  values  of  e-beam  current  density  required  to  realize  the  power  densities  can  be  estimated 
from  the  Fig.2.  It  is  evidently  that  magnitudes  of  Bmm  presented  in  the  Table  2  are  technically 
achieved.  If  the  optimal  power  spent  on  flow  ionization  presented  in  the  Table  2  will  be  hard 
to  achieve,  the  alternative  regime  with  qi<qopt ,  as  it  follows  from  the  Fig.8,  can  be  realized  at 

To  analyze  potentiality  of  MHD  control  for  improving  characteristics  of  complex  gas- 
dynamic  system  it  is  necessary  to  develop  mathematical  model  of  MHD  systems.  In  this 
chapter  we  develop  quasi- ID  model  of  MHD  generator  with  nonequilibrium  conductivity. 
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so-called  relative  power  excess  or  the  real  power  produced  by  MHD  generator  with 
nonequilibrium  conductivity.  One  can  see  that  at  some  magnitude  of  MHD  interaction 
parameter  Sv  the  //„  and  r/ion  values  become  equal.  So  self-sustained  operational  mode  for 
MHD  generator  with  nonequilibrium  conductivity  is  limited  not  only  by  parameters  of  ionizer 
and  magnetic  system,  but  also  by  the  length  of  the  MHD  generator.  Points  in  the  Fig.  10  mark 
the  limiting  values  of  MHD  interaction  parameter.  As  it  follows  from  Fig.  10  the  limiting 
MHD  interaction  parameter  for  Hall  MHD  generator  is  less  than  for  Faraday  one. 

Fig.  11  shows  relative  power  excess  produced  by  Faraday  MHD  generator  in  the 
regimes  with  £  =  |(^)(a)  and  %=0  (b).  One  can  see  that  regime  £=0  is  more  preferable 

regime  for  power  production.  While  increasing  the  Mach  number  at  MHD  generator  entrance 
the  power  excess  increases  too  in  both  regimes.  As  it  follows  from  Fig.  1  la  the  ionizer  with  a 
constant  power  density  spent  on  ionization  provide  more  power  excess  than  ionizer  with  a 
constant  e-beam  current  density  in  MHD  generator  with  £  =  For  MHD  generator  with 

£= 0  the  tendency  is  contrary.  In  comparing  Fig.  11a  and  Fig.  12  one  can  conclude  that 
characteristics  of  ionizer  and  magnetic  system,  which  according  to  (6, 12b, 24)  are 
accumulated  in  the  ratio  qi/qcr,  are  essentially  influence  on  the  power  excess  produced  by 
MHD  generator  with  nonequilibrium  conductivity. 

Thus  in  this  chapter  the  model  of  nonequilibrium  plasma  sustained  by  e-beam  is 
developed.  Simple  relations  to  calculations  of  electron  concentration  are  obtained.  Quasi- ID 
model  of  MHD  generator  with  nonequilibrium  conductivity  is  developed.  Brief  analysis  of 
the  MHD  generator  is  made.  The  main  purpose  of  MHD  generator  investigating  in  the  project 
is  control  of  flow.  Probably  maximal  power  production  is  not  necessity  for  the  purpose.  To 
formulate  requirements  to  parameters  of  the  MHD  generator  providing  optimal  control  of 
gas-dynamic  system  it  is  necessary  to  develop  mathematical  model  of  the  system.  This  will 
be  done  in  the  next  chapter. 


Fig.l.  Geometry  of  scramjet  which  is  characterized  by  the  total  turning  angle  dN=  15°  and 
Mach  designed  number  M</=10.  Characteristic  cross-sections  used  in  calculations  are  denoted 
on  the  figure. 
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3.  QUASI- 1-D  MODEL  OF  SCRAM  JET  WITH  MHD  CONTROL 


This  chapter  is  devoted  to  development  of  quasi- 1-D  mathematical  model  of 
combined  gas-dynamical  system,  which  includes  MHD  generator  with  non-equilibrium 
conductivity  and  to  determination  of  conditions  at  which  MHD  interaction  has  essential 
influence  on  performance  of  the  gas-dynamical  system.  As  the  combined  gas-dynamical 
system  let’s  analyze  the  scheme  of  scramjet  with  MHD  control  under  “AJAX”  concept  [11]. 
Simplified  scheme  of  the  scramjet  is  presented  in  the  Fig.  13.  MHD  generators,  disposed 
upstream  the  combustion  chamber,  transform  part  of  flow  enthalpy  to  electric  power.  The 
electric  power  produced  by  MHD  generators  is  transferred  to  MHD  accelerator  disposed 
downstream  the  combustion  chamber.  Part  of  the  electric  power  need  to  be  used  to  provide 
necessary  flow  conductivity  in  MHD  channel.  Part  of  produced  electric  power  can  be 
transferred  to  onboard  systems.  Influence  of  external  MHD  generator  on  flow  field  in 
scramjet  inlet  will  be  investigated  in  the  next  chapter. 

In  this  chapter  only  internal  MHD  systems  will  be  considered.  Scheme  of  scramjet 
and  some  peculiarities  of  quasi  ID  model  developed  in  the  chapter  are  shown  in  the  Fig.  14. 
We  analyze  scheme  in  which  whole  of  power  excess  produced  by  MHD  generator  with 
nonequilibrium  conductivity  is  transferred  to  MHD  accelerator.  In  the  case  when  the  MHD 
generator  will  be  used  as  a  power  source  but  not  a  means  of  scramjet  control  (not  analyzed  in 
the  project)  it  is  necessary  to  estimate  influence  of  power  extraction  (not  recovered  in  MHD 
accelerator  in  this  case)  on  the  scramjet  specific  impulse.  Flow-field  upstream  the  MHD 
generator  entrance  is  calculated  in  two-dimensional  Euler  approach.  At  cross-section  1 
(Fig.  14)  the  flow  parameters  are  averaged.  Downstream  the  cross-section  1  the  ID  approach 
is  used.  To  determine  characteristics  of  nonequilibrium  plasma  in  MHD  channel  we  use 
results  obtained  in  the  previous  chapter.  We  are  taking  into  account  that  part  r  of  the  energy 
produced  by  MHD  generator  Wion  =rWg  is  transferred  to  ionizer.  The  rest  of  the  energy, 

namely  (l -r)W  ,  is  transferred  to  MHD  accelerator.  Combustion  chamber  is  considered  in 

two  regimes  with  a  constant  pressure  along  the  channel  and  with  a  constant  gas  density  along 
the  channel.  A  mass  flow  rate  of  fuel  is  usually  much  less  than  air  mass  flow  rate  rh  ,  thus  we 
will  regard  a  fuel  supply  into  the  combustion  chamber  as  a  heat  release  without  injection  of 
mass.  Flow  in  nozzle  is  considered  as  isentropic.  Pressure  at  nozzle  outlet  is  supposed  to  be 
equal  the  static  pressure  in  incident  flow. 

We  will  use  subscripts  to  note  parameters  of  a  flow  in  corresponding  cross-sections  of 
the  engine  channel.  Let's  use  following  numbering:  0  -  incident  flow;  1  -  entrance  of  MHD 
generator;  2  -  entrance  of  the  combustion  chamber;  3  -  entrance  of  MHD  accelerator;  4  - 
entrance  of  nozzle;  5  -  exit  of  nozzle.  For  obviousness  these  designations  are  indicated  in  the 
Fig.  14.  Now  we  introduce  mathematical  models  for  subsystems  of  the  engine. 

0-1. Inlet 

The  inlet  implements  a  multi-shock  gas-dynamical  compression  of  an  incident  flow. 
The  following  characteristics  are  used:  N  -  number  of  shocks  in  an  external  part,  0^  -  a  total 
turning  angle  of  a  flow  in  the  inlet,  oin  -  a  total  pressure  recovery  coefficient.  Characteristics 
of  inlet  were  calculated  in  a  wide  range  of  inlet  configuration  for  various  values  of  flight 
Mach  number.  Typical  flow  field  in  inlet  with  N=2  in  the  case  when  flight  Mach  number  M0 
is  equal  to  design  Mach  number  Md  is  shown  in  the  Fig.  15.  Pressure  distribution  across  the 


23 


channel  for  various  flight  Mach  numbers  are  shown  in  the  Fig.  16.  Dashed  lines  in  the  figure 
are  averaged  pressures  for  corresponding  conditions. 

If  a  temperature  at  an  exit  of  the  inlet  (entrance  of  MHD  generator)  is  assigned  a  value 
Tj,  a  pressure  and  velocity  at  the  exit  are  determined  by  the  following  ratios: 

r 


(  t\ 


=  cr„ 


V^oy 


r- i 


Pi_ 

Po 

— +  c  Tj  =  — L  +  c  7J 
2  p  0  2  p  1 


(25) 


where  cp  is  specific  heat  of  air.  Total  pressure  recovery  coefficient  ain  is  calculated  from 
upper  relation  (25)  in  terms  of  calculated  numerically  ratios  pi/po  and  Tj/To,  for  any  given 
inlet. 

Fig.  17  shows  dependencies  both  of  relative  pressure  increase  pi/po  and  relative 
temperature  increase  Tj/To  in  the  inlet  with  total  turning  angle  6n=  15°  and  Mach  designed 
number  Mj=  10  upon  relative  values  of  the  inlet  throat  Fth  for  various  values  of  flight  Mach 
number. 


1-2.  MHD  generator  with  nonequilibrium  conductivity 

The  MHD  generator  is  described  in  the  approach  developed  in  previous  chapter.  It  is 
characterized  by  parameters  E,i,  ki,  enthalpy  extraction  ratio  r\  and  relative  value  of  power 
spent  on  flow  ionization  r=Wion/Wg.  The  change  of  flow  parameters  in  a  channel  of  MHD 
generator  is  determined  as  follows: 
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- 1 ,  Mach  number  at  MHD  generator  entrance  Mi  is  determined  in 
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,  factor  y/  is 


determined  as  in  (18)  but  instead  of  k  the  load  factor  ki  is  substituted. 

The  electrical  power  Wg,  produced  by  MHD  generator  is  determined  by  the  following 
ratio:  W  =  me pT0  ■  co(M 0)  •  77 .  Factor  r=Wion/Wg  which  defines  the  relative  power  spent  on 

flow  ionization  is  determined  by  the  formulae  (23).  The  power  spent  on  flow  ionization, 
ultimately,  as  a  result  of  recombination  processes,  passes  into  a  heat.  We  suppose 
approximately  that  this  additional  heat  release  is  implemented  in  the  combustion  chamber. 


2-3.  Combustion  chamber 

Set  of  equations  for  ID  model  of  combustion  chamber  in  approach  with  a  constant 
specific  heat  is  presented  here: 
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(35b) 


r+1M  >  T] 

Mf-\  AT 

One  can  see  that  inequality  (35a)  prescribes  more  intensive  increase  of  pressure  in  MHD 
generator  channel  than  inequality  (34a)  to  improve  scramjet  performance.  Fig.  19 
demonstrates  comparison  of  Maximal  Mach  numbers  for  scramjet  with  6m=\  0°  for  two 
regimes  of  combustion  chamber.  MHD  generator  with  constant  cross-sectional  area  is 
considered.  One  can  see  that  MHD  control  can  be  used  in  wider  ranges  of  flight  Mach 
numbers  for  scramjet  with  combustion  chamber  working  in  a  mode  with  a  constant  pressure. 
In  analysis  of  scramjet  the  combustion  chamber  working  in  a  mode  with  a  constant  pressure 
is  traditionally  considered.  So  hereinafter  we  will  consider  the  same  regime  in  investigating 
the  scramjet  with  MHD  control. 

To  choose  preferential  regime  of  flow  in  MHD  generator  we  compare  specific 
impulses  of  scramjet  with  MHD  control  for  various  values  of  factor.  Relative  values  of 
specific  impulse  for  scramjet  with  MHD  control  (namely  specific  impulse  of  “MHD 
scramjet”  divided  by  specific  impulse  of  classical  scramjet)  as  a  function  of  the  enthalpy 
extraction  ratio  are  shown  in  the  Fig.20a.  One  can  see  that  increase  of  factor  ^  causes  the 
scramjet  specific  impulse  to  be  increased.  Fig.20b  demonstrates  dependency  of  relative 
power  spent  on  ionization  r=Wi0n/Wg  as  a  function  the  enthalpy  extraction  ratio.  While 
increasing  the  %i  the  factor  r  rise  take  place.  The  flow  regime  characterized  by  minimal  £/  is 
preferable  if  MHD  generator  with  nonequilibrium  conductivity  is  considered  as  a  power 
source.  Requirements  for  ^  parameter  in  the  case  when  MHD  generator  is  used  for  scramjet 
control  are  opposite.  In  this  research  we  consider  possibilities  of  nonequilibrium  MHD 

generator  for  flow  control  and  so  maximal  realizable  factor  £/  is  required.  Factors  <^j  >  ( i//, ) 
correspond  to  convergent  channel  and  undesirable  to  avoid  additional  shock  waves.  So 
hereinafter  we  will  consider  MHD  flows  characterized  by  factor  which 

approximately  describes  MHD  generator  with  constant  cross-sectional  area. 

In  the  calculations  the  factor  was  assumed  to  be  zero.  Numerical  calculations  of 
specific  impulse  for  various  values  of  factor  have  shown  that  flow  regime  in  channel  of 
MHD  accelerator  practically  doesn’t  influence  on  the  specific  impulses  value.  So  hereinafter 
we  suppose  £j=0. 

Various  regimes  of  power  spending  on  ionization  of  flow  are  compared  in  the  Fig.21. 
One  can  see  that  ionization  regime  with  a  constant  power  density  spent  on  ionization  gives  us 
opportunity  to  reach  greater  values  of  specific  impulse  and  power  excess  than  the  regime  with 
a  constant  e-beam  current  density.  This  tendency  is  correct  both  for  Faraday  and  Hall  MHD 
generators.  So  hereinafter  we  will  consider  the  regime  with  a  constant  power  density  spent  on 
ionization. 

Fig.22  demonstrates  dependencies  of  relative  specific  impulse  as  a  function  of  the 
load  factor  for  various  flight  Mach  numbers  and  various  values  of  enthalpy  extraction  ratio. 
One  can  see  that  positive  influence  of  MHD  interaction  on  specific  impulse  increase  while 
decreasing  the  flight  Mach  number. 

Fig.23  demonstrates  dependencies  of  MPCE  specific  impulse  upon  the  MHD 
interaction  parameter  for  various  relative  values  of  inlet  throat.  Results  are  presented  for 
limiting  case  for  which  qi/qcr  =0.  One  can  see,  the  more  is  the  inlet  throat  Fth  the  greater  is  the 
positive  influence  of  MHD  interaction  on  specific  impulse  of  scramjet  both  for  Faraday  and 
Hall  MHD  generators.  Influence  of  MHD  generator  type  on  specific  impulse  of  scramjet  is 
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Fig.  16.  Relative  pressure  (a)  and  temperature  (b)  distributions  in  scramjet  inlet  at  cross-section  (1) 
for  various  flight  Mach  numbers  (shown  near  the  curves)  obtained  in  2D  calculations  (solid  lines) 
and  used  in  ID  model  (dashed  lines).  Md=10,  0n=15°,  Fth=0.12. 
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Fig.27  Relative  power  excess  produced  by  MHD  generator  in  scramjet  vs  MHD  interaction 
parameter  for  various  inlet  configurations  and  various  load  factors.  M,/=10  F,/;=0.1,  Mo=6, 
q/qcr= 0.05,  ks=  1.2.  (a)  Faraday  MHD  generator,  (b)  Hall  MHD  generator  with  p=2 
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Fig.29b  Relative  specific  impulse  of  scramjet  with  MHD  control  as  a  function  of  MHD  interaction 
parameter  for  various  load  factor  of  MHD  accelerator  ks\  a)  Faraday  MHD  generator,  b)  Hall  MHD 
generator  with  p=2.  Mcj=  1 0  F,/,=0.12,  Mo=6,  q,/qcr=0.05. 
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Functions  fx,  fy,  and  qg  from  the  right  part  of  (36)  are  determined  from  (37)  as  a 
functions  of  parameters  Bx,  By,  vx,  vy,  Ex,  Ey,  Ez  and  cr.  We  consider  ideally  sectioned  Faraday 
MHD  generator  for  which  current  density  has  only  z  component  (j  =  ( 0, 0,  j_ ) ).  In  this  case 
MHD  interaction  can  be  characterized  by  dimensionless  factor  k=-Ez/voB0 by  the  ratios: 


Computational  procedures  for  solution  the  set  of  equations  (36)  are  based  on  the 
explicit  high-resolution  shock-capturing  Godunov-type  scheme  for  space-marching 
calculations  of  stationary  supersonic  flows.  Initial  conditions  are  given  by  incoming  flow. 
Normal  velocity  component  on  the  wall  is  assumed  to  be  zero. 

Fig.35  demonstrates  pressure  contours  in  scramjet  inlet  without  MHD  interaction  at 
design  conditions.  Fragment  of  computational  grid  from  red  rectangular  area  marked  in  the 
Fig.35a  are  shown  in  the  Fig.35b.  In  this  calculation  the  size  of  grid  was  2529x200.  Grid  is 
adaptive  in  x  direction.  Fragments  of  computational  grid  in  MHD  controlled  inlet  are  shown 
in  the  Fig.36.  Fig.36a  corresponds  to  inlet  without  MHD  interaction.  In  this  case  the  grid  is 
uniform.  Fig.36b  corresponds  to  inlet  with  MHD  generator  located  in  the  range  I  <x<I .2.  In 
this  case  the  grid  becomes  nonuniform.  Step  of  the  grid  decreases  in  x  direction  in  the  MHD 
interaction  area. 


In  calculations  of  MHD  controlled  inlet  we  will  consider  various  simple  geometries  of 
ionized  region,  which  are  shown  in  the  Fig.37.  MHD  generator  with  nonequilibrium 
conductivity  is  considered.  Flow  conductivity  a  in  ionized  region  is  a  function  of  flow 
parameters  and  power  density  spent  on  flow  ionization  qt  (the  chapter  2).  Gas-dynamic  flow 
compression  in  scramjet  inlet  occurs  in  series  of  oblique  shocks.  So  flow  parameters  in 
scramjet  without  MHD  interaction  (see  Fig.35a)  are  inhomogeneous.  Thus  requirements  to 
parameters  of  ionizer  and  magnetic  system,  which  provide  self-sustained  operational  mode 
for  MHD  generator  with  nonequilibrium  conductivity  will  be  a  function  of  geometry  of  inlet 
and  position  of,  ionized  region.  In  the  general  case  the  condition  for  self-sustained 
operational  mode  for  the  MHD  generator  can  be  written  in  the  form: 

j] \qi{x,y)dxdy  g(qi,x,y)dxdy  (39) 

V  V 


where  V  is  a  volume  in  which  ionization  of  a  flow  and  extraction  of  electric  power  by  MHD 
generator  take  place,  qg  is  calculated  according  to  (38).  In  calculations  we  suppose  that  qt  has 
a  constant  value  in  the  region  of  ionization.  Function  qg  depends  on  the  x  and  y  coordinates 
because  variables  vx,  vy  and  cr  are  functions  of  the  coordinates.  The  magnitude  of  qt  at  which 
in  expression  (39)  is  implemented  a  sign  of  equality,  corresponds  to  critical  operational  mode 
of  MHD  generator  with  nonequilibrium  conductivity.  It  is  easy  to  see  from  (39)  that  the 
following  equation  determines  critical  power  in  nonlocal  approach: 

qcr  =]r\\(lg  (Qcr  >  T  y)dxdy 

*  V 
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coordinates  of  the  upper  border  of  the  region  Y2.  Conditions  are  the  same  as  in  the  Fig.61. 
Similarly  to  the  air  mass  flow  in  inlet  the  flow  compression  increases  while  increasing  the  Y2 
value.  Positions  of  bottom  border  Yj  for  the  ionized  region  optimal  for  increasing  the  flow 
compression  are  differ  from  positions  optimal  for  increasing  the  air  mass  flow.  According  to 
the  Fig.62,  values  of  7;,  optimal  for  increasing  of  flow  compression,  are  less  than  0.7m. 
Fig.63  demonstrates  how  cp/cpo  and  n/no  depend  on  the  orientation  of  magnetic  induction  for 
various  position  of  ionized  region  in  x-direction.  The  ionized  region  has  rectangular  shape 
and  it  is  located  in  Xi<x<Xi+I  and  0.9<y<1.0.  One  can  see  that  for  this  configuration 
maximal  air  mass  flow  and  flow  compression  are  achieved  at  a«105°.  Air  mass  flow  in  the 
inlet  not  essentially  depends  on  the  X )  value.  On  the  contrary  the  flow  compression  in  the 
inlet  significantly  depends  on  the  Xj.  So  we  can  conclude  that  shape  and  position  of  ionized 
region  are  very  important  factors  which  influence  on  characteristics  of  MHD  controlled  inlet. 

And  now  we  consider  influence  of  B  and  q,  values  on  characteristics  of  MHD 
controlled  inlet  in  the  case  when  ionized  region  is  the  band  located  in  3<x,<4.  Fig.64  shows 
the  ratio  (p/cpo  as  a  function  of  magnetic  induction  B  for  various  values  of  qt.  It  is  easy  to  see 
that  there  are  optimal  values  of  magnetic  induction,  depending  on  the  power  density  spent  on 
flow  ionization,  at  which  air  mass  flow  in  MHD  controlled  inlet  achieve  maximum.  The 
optimal  magnetic  induction  decreases  while  increasing  the  qt.  Maximal  normalized  relative 
air  mass  flow  decreases  while  increasing  the  power  density  spent  on  ionization  in  the  range 
from  O.lW/cm3  to  10W/cm3.  So  to  achieve  maximal  air  capture  in  MHD  controlled  inlet  it  is 
necessary  to  ensure  ionization  with  q,  not  greater  than  1  W/cm3 .  Fig.65  demonstrates  the  ratio 
7t/7to  as  a  function  of  magnetic  induction  B  for  various  values  of  q,.  Similarly  to  air  mass  flow 
the  flow  compression  achieves  a  maximal  value  at  some  optimal  magnetic  induction  and  the 
optimal  value  of  magnetic  induction  is  decreasing  function  of  qt  value.  The  maximal  flow 
compression  in  the  inlet  is  achieved  at  maximal  power  density  spent  on  ionization.  Thus  both 
the  power  density  spent  on  ionization  and  the  value  of  magnetic  induction  are  important 
parameters,  which  determine  characteristics  of  MHD  controlled  inlet. 

Obtained  results  show  that  MHD  control  allows  us  to  modify  flow  field  in  inlet  of 
scramjet  at  off-design  conditions.  By  using  MHD  influence  on  flow  we  can  control  oblique 
shocks  position,  air  mass  flow  and  flow  compression  in  inlet.  In  proper  choice  of  region  of 
MHD  interaction,  configuration  of  magnetic  field  and  load  factor  the  MHD  control  in 
scramjet  inlet  allows  us  essentially,  improve  inlet  performance  at  off-design  conditions. 


Fig.34  Geometry  of  MHD-controlled  inlet  with  arbitrary  region  of  flow  ionized  by  e-beam, 

$v=15°,  Mq=l0,  Fth=0.l. 
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c) 

Fig.39.  Pressure  colors  maps  in  stream  flowing  the  wedge  with  0=10°.  Mo=8,  ^,=4W/cm  , 
A=0.5.  a)  B=0,  b)  B=3T,  a=45°,  c)  B=3T,  a=90°.  Ionized  region  located  in  2<x<2.5. 
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1.1 


a) 


Fig.41  Temperature  contours  in  stream  flowing  the  wedge  with  0=10°.  M0=8,  c/,=4W/cm  , 
A=0.5.  a)  B=3T,  a=45°,  b)  B=3T,  a=90°.  Ionized  region  located  in  2<x<2.5. 
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Fig.44  Temperature  contours  in  stream  flowing  the  wedge  with  0=10°.  Mo=8,  </,=4W/cm3, 
£=0.5.  B=3T,  a=135°.  Ionized  region  located  in  2<x<2.5. 


Fig.45.  Position  of  oblique  shock  in  stream  flowing  around  the  wedge  with  0=10°.  Mo=8, 
<^=4W/cm3,  £=0.5.  Red  contours  correspond  to  wedge  without  MHD  control;  for  blue 
contours  B=5T,  a=135°;  for  green  contours  B=1T,  a=45°.  Ionized  region  located  in  2<x<2.5. 
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Fig.48  Density  distribution  across  the  channel  at  x=7. 1  in  MHD  controlled  inlet  at  various 
values  of  power  density  spent  on  ionization.  M(j=\  0,  0N=\5°,  Fth= 0.12,  B=  IT,  a=45°,  k=0.5, 
qi  values  are  shown  in  the  figure.  Ionized  region  is  located  in  0<x,<5,  y<0.2x. 


Fig.49.  Mach  number  contours  in  scramjet  inlet.  Mo= 8,  Mj=  10,  8n=\  5°,  F, /,=().  1 2. 


66 


q. ,  W/cm3 


b)  q, ,  W/cm3 

Fig.51.  Relative  mass  flow  rate  (a)  and  pressure  increase  (b)  in  MHD  controlled  inlet  vs 
power  density  spent  on  ionization.  Ionized  region  is  located  in  xi<x,<xi+l.  Mo=8,  Md=  10, 
0n=  15°,  F^=0.12,  B=  IT,  &=0.5,  a =135°,  x;  values  are  shown  in  the  figure. 
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Fig.62.  Normalized  relative  pressure  increase  in  MHD  controlled  inlet  vs  7;  for  various 
values  of  Y2.  Mo= 8,  Mj=  10,  $v=15°,F,/,=0.12,  B= 3T,  a=135°,  ^=0.1W/cm3,  k=0.5.  Ionized 
region  is  located  in  3<x,<4  and  Yi<y<Y2 
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Fig.63.  Normalized  relative  air  mass  flow  and  pressure  increase  in  MHD  controlled  inlet  vs  a 
for  various  values  of  20.  Mo= 8,  M,/=10,  0N=l5°,Fth=OA2,  B= 3T,  q,= lW/cm3,  k=0.5.  Ionized 
region  is  located  mXi<x<Xi+l  and  0.9<y<1.0. 
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B,  T 

Fig.64.  Normalized  relative  air  mass  flow  in  MHD  controlled  inlet  vs  magnetic  induction  B 
for  various  values  of  qt.  Mo= 8,  M(/=10,  #v=15°,F,/,=0.12,  a=135°,  k=0.5.  Ionized  region  is 
located  in  3<x<4. 


Fig.65  Normalized  relative  pressure  increase  in  MHD  controlled  inlet  vs  magnetic  induction 
B  for  various  values  of  qt.  Mo= 8,  Md=  10,  6V=15°,F(/,=0.12,  a=135°,  k=0.5.  Ionized  region  is 
located  in  3<x<4. 
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5.  CONCLUSIONS 


The  main  results  obtained  in  this  report  are  as  follows: 

•  The  model  of  nonequilibrium  plasma  sustained  by  e-beam  is  developed.  Simple 
relations  to  calculations  of  electron  concentration  as  a  function  of  power  density  spent 
on  ionization  and  flow  parameters  are  obtained. 

•  Quasi- ID  model  of  MHD  generator  with  nonequilibrium  conductivity  is  developed. 
Requirements  to  parameters  of  ionizer  and  magnetic  systems  at  which  nonequilibrium 
MHD  generator  operates  in  self-sustained  mode  are  formulated. 

•  Quasi- ID  model  of  scramjet  with  MHD  control  by  internal  flow  is  developed. 
Requirements  for  parameters  of  MHD  generator  and  other  subsystems  of  the  scramjet, 
at  which  MHD  control  increases  scramjet  specific  impulse,  are  formulated. 
Limitations  on  flight  Mach  number  at  which  internal  MHD  generator  can  improve 
scramjet  performance  are  obtained.  Analysis  of  scramjet  with  MHD  control  shows 
that  MHD  generator  with  nonequilibrium  conductivity  can  be  used  on  hypersonic 
aircraft  both  for  electric  power  production  and  for  scramjet  specific  impulse 
increasing.  Requirements  for  parameters  of  the  MHD  generator  depend  on  its 
functionality. 

•  2-D  model  of  MHD  controlled  inlet  of  scramjet  is  developed.  Numerical  calculations 
of  flow  fields  in  MHD  controlled  inlet  are  carried  out  in  a  wide  range  of  parameters 
variation.  Obtained  results  show  that  MHD  control  allows  us  to  modify  flow  field  in 
inlet  of  scramjet  at  off-design  conditions.  By  using  MHD  influence  on  flow  we  can 
control  oblique  shocks  position,  air  mass  flow  and  flow  compression  in  inlet.  In 
proper  choice  of  region  of  MHD  interaction,  configuration  of  magnetic  field  and  load 
factor  the  MHD  control  in  scramjet  inlet  allows  us  essentially,  improve  inlet 
performance  at  off-design  conditions. 

Next  steps  need  to  be  done  in  order  to  create  more  overall  physical  model  of  MHD 
controlled  inlet.  It  is  necessary  to  keep  in  mind  real  configurations  of  magnetic  field,  consider 
real  propagation  of  e-beam  in  electrical  and  magnetic  fields  in  nonhomogeneous  medium 
taking  into  account  processes  of  plasma  kinetics.  Influence  of  space  charge  distribution  on 
electric  field  and  e-beam  propagation  need  to  be  considered  too. 
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DESCRIPTION  OF  VARIABLES 

final  report  on  the  ISTC  project  N°  2088p,  “Theoretical  and  experimental  research  of 

capabilities  of  MHD  technology  to  control  gas  flow  with  non-equilibrium  ionization” 

A 

Cross-sectional  area 

2 

m 

a, ,  bi 

Fitting  parameters 

B 

Magnetic  induction 

T 

Bcr 

Critical  value  of  magnetic  induction 

T 

Bmin 

Minimal  value  of  magnetic  induction 

T 

cp 

Specific  heat  at  constant  pressure 

J/(kgK) 

cv 

Specific  heat  at  constant  volume 

J/(kgK) 

E 

Electric  field  intensity 

V/m 

e 

Electron  charge 

C 

Eb 

Initial  energy  of  electrons  in  e-beam 

keV 

Ff 

Lorentz  force 

N/m3 

Fth 

Relative  value  of  the  inlet  throat 

_g_ 

Acceleration  of  gravity 

m/s2 

Sen 

Averaged  relative  velocity  of  electrons 

m/s 

h 

Static  enthalpy 

J/kg 

Hu 

Calorific  value  of  fuel 

Mg 

I 

E-beam  induced  ionization  rate 

-3  -1 

cm  s 

Ep 

Specific  impulse 

s 

j 

Current  density 

A/m2 

jb 

E-beam  current  density  A/cm2,  A/m 2 

k 

Load  factor 

ka 

Rate  constant  for  attachment  of  electrons 

cm6/s 

kc 

Electron  scattering  rate  constant 

m3/s 

kd 

Rate  constant  for  detachment  of  electrons 

cm3 Vs 

L 

Length  of  MHD  generator 

m 

L0 

Stochiometric  factor 

M 

Mach  number 

M  ,M0 

Flight  Mach  number 

Md 

Design  Mach  number 

Mmax 

Maximal  flight  Mach  number.  (MHD  interaction  increase  specific 

me 

impulse  of  scramjet  in  the  case  when  M  <Mmax) 

Electron  mass 

kg 

m 

Air  mass  flow 

kg/s 

N,  n 

Gas  concentration 

-3 

cm 

n. 

Negative  ions  concentration 

-3 

cm 

n+ 

Positive  ions  concentration 

cm 3 

ne 

Electron  concentration 

-3 

cm 

No2 

Concentration  of  oxygen  molecules  in  air 

-3 

cm 

P 

Pressure 

Pa,  atm 

Qcr 

Critical  power  density  spent  on  flow  ionization 

W/cm3, 

q, g 

Power  density  produced  by  MHD  generator 

W/m3 

W/cm3, 

Vi 

Power  losses  by  electron  beam  in  unit  of  time  in  unit  of  volume 

W/m3 

W/cm3, 

(power  density  spent  on  ionization) 

W/m3 

Qopt 

Optimal  value  of  the  power  density  spent  on  flow  ionization 

W/cm3 

qr 

Power  density  released  in  result  of  plasma  recombination 

W/m3 

Qen 


R 

Sv 

T 

t 

Te 

Tg 

V 

V 
Wo 
Wa 

Wg 

Wi 

wion 

x,y,z 

V 


ii 


>  S 
ion 


0 

N 


0 

N 


in 


Averaged  transport  cross-section  of  electron  collisions 


m 


2 


Ratio  of  power  spent  on  ionization  of  flow  in  MHD  generator  and 
power  produced  by  MHD  generator  (r=  Wion/Wg) 

Gas  constant 

MHD  interaction  parameter 

Temperature 

Time 

Electron  temperature 
Gas  temperature 
Flow  velocity 
Volume 
Energy  of  flow 

Power  transferred  to  MHD  accelerator 
Power  produced  by  MHD  generator 
Ionization  cost 

Power  put  to  flow  ionization  in  channel  of  MHD  generator 
Cartesian  coordinates 
Electron  stopping  power  in  air 
Air-fuel  ratio  (used  in  chapter  3) 

Angle  of  magnetic  field  orientation  (used  in  chapter  4) 

Rate  constant  for  electron-ion  recombination  (used  above  equation  14) 
Hall  parameter  (used  after  equation  14  ) 

Rate  constant  for  ion-ion  recombination 
Specific  heat  ratio 
Enthalpy  extraction  ratio 

Relative  power  spent  on  flow  ionization  (  ion=Wi0n/Wo) 

Relative  air  mass  flow  in  scramjet  with  MHD  control 
Relative  air  mass  flow  in  scramjet  without  MHD  control 
Nozzle  non-ideality  factor 
Electron  mobility 

Flow  compression  in  inlet  with  MHD  control 
Flow  compression  in  inlet  without  MHD  control 
Total  turning  angle  of  flow  in  scramjet  inlet 

Gas  density 

Electrical  conductivity 

Total  pressure  recovery  coefficient 

Factor  which  determines  regime  of  MHD  flow 

Factor  which  characterizes  the  MHD  generator  type  (Faraday  or  Hall) 


J/(kgK) 

K 

s 

K 

K 

m/s 

3 

m 

W 

W 

W 

eV 

W 

Me  V  cm2 g1 

degree  of 
arc 
an  /s 

cm3/s 


m2/(V  s) 


degree  of 
arc 

g/cm3,  kg/m3 
mho/m 


Subscripts 

x,  y,  z  Denotes  the  vector  projection  onto  corresponding  axes  of 

Cartesian  coordinates 

0,1, 2, 3, 4,5  Denotes  the  parameters  at  corresponding  cross-sections  of 

scramjet  with  MHD  control  according  to  the  Fig.  14. 
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which  self-sustained  operational  mode  of  ionizer  and  MHD  generator  is  realized  will  be  formulated.  Possibilities  of  using  of  MHD  control  in 
gas-dynamical  systems  will  be  considered.  Traditional  use  of  MHD  transformation  of  energy  to  produce  electric  power  is  well  understood  both 
theoretically  and  experimentally  nowadays.  Nontraditional  use  of  MHD  to  control  flow  parameters  and  hence  performance  of  gas-dynamical 
systems  that  include  MHD  system  as  a  subsystem  is  not  well  understood.  MHD  interaction  acts  as  bulk  force  and  power  upon  flow  and  allows 
varying  flow  parameters.  In  this  case  no  mechanical  parts  are  used  to  control  the  flow  parameters.  Therefore,  this  method  of  control  can  be 
considered  as  a  perspective  one  for  gas-dynamical  system  with  a  fixed  geometry  especially  for  high  enthalpy  flows.  As  a  possible  range  of 
application  for  MHD  control,  two  types  of  gas-dynamical  system  can  be  considered:  a)  gas-dynamical  systems  which  must  ensure  variable 
parameters  of  flow  at  outlet  of  system  with  fixed  parameters  at  inlet  of  one  and  b)  gas-dynamical  systems  which  must  ensure  limited  range  of 
variation  for  some  flow  parameters  at  outlet  of  system  with  variable  parameters  at  inlet.  The  main  requirement  to  realize  MHD  influence  on 
flow  is  flow  conductivity.  In  most  gas-dynamical  systems  thermal  ionization  of  flow  is  inappreciable  and  therefore  it  is  necessary  to  use 
methods  of  non-equilibrium  plasma  creation.  Various  methods  of  flow  ionization  will  be  analyzed  in  the  project.  To  ensure  autonomous 
regime  of  system  of  MHD  control  self-sustained  operational  mode  of  ionizer  and  MHD  generator  is  considered.  Power  spent  on  flow 
ionization  is  expected  to  be  less  than  power  produced  by  MHD  generator  for  this  operational  mode. 

In  process  of  implementation  of  the  project  the  next  activities  will  be  carried  out:  theoretical  research  of  methods  of  non-equilibrium  plasma 
creation;  development  of  mathematical  model  of  an  ionizer;  numerical  calculations  of  electrical  conductivity  of  plasma;  determination  of  range 
of  parameters,  at  which  self-sustained  operational  mode  of  an  ionizer  and  MHD  generator  exists;  development  of  mathematical  models  of 
gas-dynamical  systems  which  include  ionizer  and  MHD  generator  as  subsystems  in  quasi-1-D  and  2-D  approaches;  numerical  calculations  of 
flow-fields  in  the  gas-dynamical  systems;  determination  of  conditions  at  which  MHD  control  appreciably  modifies  characteristics  of 
investigated  systems;  development  and  creation  of  experimental  model  of  ionizer;  experimental  testing  and  modification  of  mathematical 
model  of  ionizer;  theoretical  research  of  influence  of  configuration  of  magnetic  field  and  various  types  of  MHD  generator  on  flow-field  and 
performance  of  gas-dynamical  system;  formulating  recommendations  to  the  choice  of  parameters  of  ionizer  and  MHD  generator  to  provide 
maximal  influence  on  flow-field  in  gas-dynamical  system;  development,  creation  and  test  of  model  which  includes  MHD  system  for  flow 
control. 


a  form: 


Set  of  equations  for  stationary  MHD  flow  in  quasi- ID  approach,  according  to  [7]  has 
/?vd=const 


dv  dp  -  - 

pv—+—=fx  =  (j*b)x 

dx  dx 


.2\ 


h  + 


(13) 


pv- 


=  jE 


dx 
P  =  RpT, 

where  p  is  the  flow  density,  v  is  the  flow  velocity,  p  is  the  static  pressure,  h  is  the  enthalpy,  A 
is  the  MHD  channel  cross-sectional  area,  x  is  the  longitudinal  coordinate,  B  is  the  magnetic 
field,  E  is  the  electric  field,  j  is  the  current  density,  F  is  the  Lorentz  force  We  assume  that 

B  =  (o,  5,0) ,  E  =  ( Ex , 0,  E_ )  and  v  =  (vx,0,0).  The  generalized  Ohm’s  law  determines 
relations  between  electromagnetic  components: 

(14) 


i  +  p  (y  x  B  j  =  <7  ^E  +  v  x  B  j 


J 

where  p  is  the  electron  mobility.  In  considered  configuration  of  magnetic  field  the  current 
density  is  determined  from  equation  (14)  by  the  ratio: 


J  = 


°(E.+ PE..+ PB,v,)  <T  (PE,-E,-B,v,) 

9  ^9 


i +/r 


l +/r 


where  fd  is  the  Hall  parameter 


v  J 

J3  =  p  B  .  We  will  consider  two  configurations  for  MHD  generator.  The  first  is  ideally 

sectioned  Faraday  MHD  generator,  for  which  j  =  ( 0, 0,  j_ ) ,  E  =  ( Ex ,  0,  E_ ) .  In  this  case  we 

obtain  the  following  expressions  for  the  right  part  of  equations  (13): 

(j'xB)  =  (/x5)  =-(l -k)aB2v 

^  ,  (15a) 

(y-£)  =  -A:(l-/t)cr52v2 

where  the  load  factor  k  is  determined  by  the  ratio  k  =  -E_j(vx  •  Bv ) . 

The  second  configuration  is  the  Hall  MHD  generator,  for  which  is  =  (isx,0,0)  and 
j  =  ( jx ,  0,  jz ) .  In  this  case  we  obtain: 

(7x5)^=-(1  +  ^2)(t52v/(1  +  ^2) 

(}  •  1)  =  -k  (1  -  k)  <7B2v2p2  j  (\  +  P2 ) 1 ’ 

where  the  load  factor  k  is  determined  by  the  ratio  k  =  -Ex  /( /dvxBy ) . 

In  the  Hall  MHD  generator  Lorentz  force  has  not  only  x  component  but  also  z 
component  and  ratio  of  the  components  is: 

\FJFPP(\-k)l{\  +  kp2)  (16) 

Therefore  ID  approach  for  Hall  MHD  generator  is  quite  correct  only  when  \FZ  / Fx\  <k  1 . 


(15b) 
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Fig.7.  Critical  values  of  magnetic  induction  calculated  in  different  approaches  for  various  values 
of  flight  Mach  numbers.  MHD  generator  is  located  in  cross-section  2  of  scramjet  (Fig.l).  Free 

stream  dynamic  pressure  is  equal  to  40  kPa. 

IQ-, - 


Fig.8.  Range  of  values  of  magnetic  induction  B  and  power  spent  on  flow  ionization  qt,  which 
ensure  the  self-sustained  operational  mode  of  MHD  generator  with  non-equilibrium  conductivity 
with  MHD  interaction  parameter  ^>0.1.  MHD  generator  is  located  in  cross-section  1,  Mx=6. 
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compared  in  the  Fig.24a.  For  small  value  of  Sv  the  Hall  MHD  generator  provides  more 
noticeable  increase  of  scramjet  specific  impulse  than  Faraday  one.  For  SV>0.7  the  Faraday 
MHD  generator  gives  greater  increment  of  scramjet  specific  impulse  than  Hall  one.  Fig.24b 
shows  dependencies  of  the  enthalpy  extraction  ratio  for  various  relative  values  of  inlet  throat. 
One  can  see,  the  more  is  the  inlet  throat  Fth  the  greater  is  the  value  rjg.  Faraday  MHD 
generator  produces  more  electric  power  than  Hall  one  in  all  the  cases. 

Fig.25a  demonstrates  how  relative  specific  impulse  of  scramjet  with  MHD  control 
depends  upon  MHD  interaction  parameter  for  various  relative  power  density  spent  on 
ionization  q/qcr-  One  can  see  that  increase  of  q/qcr  value  causes  the  scramjet  specific  impulse 
to  be  decreased  both  for  Faraday  and  Hall  MHD  generators.  In  the  case  of  q/qcr  >0  there  are 
optimal  MHD  interaction  parameters  at  which  relative  specific  impulse  has  a  maximum.  The 
maximal  impulse  and  optimal  MHD  interaction  parameter  decrease  while  increasing  the  ratio 
q/qcr-  The  same  behavior  is  observed  for  power  excess  presented  in  the  Fig.25b. 

Dependencies  of  relative  specific  impulse  upon  MHD  interaction  parameters  for 
scramjet  with  two  different  total  turning  angles  in  inlet  are  shown  in  the  Fig.26.  One  can  see 
that  positive  effect  of  MHD  control  in  scramjet  increase  while  decreasing  the  total  turning 
angle.  This  tendency  is  observed  for  various  values  of  the  load  factor  both  for  Faraday  MHD 
generator  (Fig.26a)  and  Hall  MHD  generator  (Fig.26b).  The  same  behavior  is  observed  for 
power  excess  presented  in  the  Fig.27.  Besides  that,  one  can  see  that  at  small  MHD  interaction 
parameter  maximal  power  generation  is  achieved  at  the  load  factor  ki= 0.5.  For  Faraday  MHD 
generator  (Fig.27a)  optimal  load  factor  becomes  greater  than  0.5  while  increasing  Sv 
parameter.  For  Hall  MHD  generator  (Fig.27b)  opposite  tendency  is  observed,  namely  the 
optimal  load  factor  becomes  less  than  0.5  while  increasing  Sv. 

Figs.28  demonstrate  how  specific  impulse  and  power  characteristics  of  scramjet  with 
Hall  MHD  generator  depend  on  the  Hall  parameter.  As  it  follows  from  the  Fig28a  the 
scramjet  specific  impulse,  at  small  values  of  Sv  parameter  monotonically  decreases  while 
increasing  the  Hall  parameter.  At  Sv  parameter  near  to  unity  the  dependence  becomes  non¬ 
monotonic.  The  relative  power  excess  (Fig.28b)  monotonically  increases  while  increasing  the 
Hall  parameter.  One  can  see  that  both  specific  impulse  and  power  excess  for  scramjet  with 
Hall  MHD  generator  are  trending  toward  corresponding  dependencies  for  scramjet  with 
Faraday  MHD  generator  while  increasing  the  Hall  parameter. 

Figs.29  demonstrate  influence  of  MHD  accelerator  load  factor  on  the  relative  specific 
impulse  of  scramjet  with  MHD  control.  One  can  see  that  decrease  of  the  load  factor  of  MHD 
accelerator  ks  causes  the  scramjet  specific  impulse  to  be  increased  both  for  scramjet  with 
Faraday  and  Hall  MHD  generators.  This  behavior  is  observed  for  various  values  of  load 
factor  of  MHD  generator  ki,  but  for  ki= 0.5  the  influence  of  k3  on  specific  impulse  is  more 
significant.  It  is  evidently  that  maximal  specific  impulse  can  be  obtained  at  ks=  1,  but  it  is 
unachievable  limit,  because  the  power  density  input  in  MHD  accelerator  qa  is  proportional  to 
k^l-ki).  Thus,  when  k3  — >  1  the  volume  and  hence  the  length  of  MHD  accelerator  La  — »  oo  . 

In  reality  it  is  inexpedient  to  take  up  the  kj  factor  less  than  1.1. 

Fig.30  demonstrates  contours  of  relative  specific  impulse  for  scramjet  with  MHD 
control  for  various  values  of  q/qcr  ratio  in  coordinates:  the  load  factor  ki  and  the  MHD 
interaction  parameter  Sv.  Left  contours  in  the  Fig.30  correspond  to  Faraday  MHD  generator 
and  right  contours  correspond  to  Hall  MHD  generator  with  p=2.  Values  of  q/qcr  ratio  are 
shown  near  the  corresponding  figures.  One  can  see  that  contours  for  Hall  MHD  generator  are 
quite  similar  to  mirroring  of  Faraday  MHD  generator  contours.  When  q/qcr  =0,  maximal 
specific  impulse  is  achieved  at  intermediate  values  of  load  factor  kj  and  maximal  values  of 
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Fig.25  Relative  specific  impulse  (a)  and  power  excess  (b)  in  scramjet  with  MHD  control  as  a 
function  of  MHD  interaction  parameter  for  various  relative  values  of  q/qcr.  6n=  15°,  F,/,=0.12, 
Md=  10,  M0= 6,  k]= 0.5,  *j=1.2. 
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a)  b) 

Fig.31  Contours  of  relative  specific  impulse  (a)  and  power  excess  (b)  for  scramjet  with  Faraday 
MHD  generator.  0N=\5°,  Md=  10,  F,/,=0.12.  Mo= 6,  L=2m,  ki=0.5,  ks=l.l 


a)  b) 

Fig.32  Contours  of  relative  specific  impulse  (a)  and  power  excess  (b)  for  scramjet  with 
Faraday  MHD  generator.  6C=1 5°,  Mj=  10,  C,/,=0. 12.  Mq=6,  L=2m,  qi=lW/cm3,  ks=l.l 
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